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This paper is concerned with the distortion produced in a frequency 
modulation wave when it passes through a filter. The phase or frequency 
modulation representing the signal is assumed to be a band of gaussian noise. 
The main result is an expression for the power spectrum W t (f) of the output 
phase angle d(t). This expression holds for any filter, contains all of the 
distortion terms due to second and third order modulation, and is suited 
to computer evaluation. It is useful in many cases, but it has the shortcoming 
of not containing any modulation terms higher than the third order. 

A second result is an approximation to W e (J), based on log (1 + x) tt x 
{that is, a " first-order" approximation), which is encountered in the deriva- 
tion of the main result. Although it does not contain all of the second and 
third order modulation terms, it does contain higher order modulation terms 
which may give most of the distortion in some cases. The results given here are 
compared with those obtained earlier. 

I. PREFACE 

This work is a sequel to "Distortion and Crosstalk of Linearly 
Filtered Angle-Modulated Signals" by E. Bedrosian of the Rand 
Corporation and myself. 1 One of the principal results of that paper 
is an expression for the distortion produced when a frequency modula- 
tion wave, modulated by gaussian noise, passes through a filter as- 
sumed to be symmetrical about the carrier frequency. 

The assumption of symmetry simplified the analysis, but led to 
zero second order modulation; and consequently the results do not 
apply to many cases of practical interest. 

The main result of this paper is an expression for the distortion 
which contains all of the second and third order modulation terms pro- 
duced by a general filter. It includes the earlier result as a special case. 
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As before, the input phase angle is assumed to be a gaussian noise. 
Some time ago Mr. Bedrosian and I worked out, independently, 
the second order modulation terms. His analysis is somewhat different 
from mine and throws a different light on the problem. Since each ap- 
proach is of interest in its own right and since we were unable to 
combine the two without losing useful results, we decided to publish 
our work separately. An early version of Bedrosian's analysis is given 
in a RAND memorandum. 2 

II. INTRODUCTION 

When an angle-modulated wave (FM or PM) passes through a 
filter, the signal becomes distorted. For a multichannel system this 
distortion may produce crosstalk. In many practical cases the second 
and third order modulation terms give a good measure of the distor- 
tion. These terms have been studied by a number of investigators. In 
this paper we obtain some general expressions for them for the case in 
which the modulation is gaussian. 

Our main results includes all of the second and third order modu- 
lation products. In this respect, it is more general than some of the 
earlier expressions for the distortion (see Medhurst, 3 Magnusson, 4 and 
Liou 5 ) . However, it does not give higher order modulation terms, some 
of which appear in earlier "first order" approximations. A first order 
approximation (similar to the earlier ones) occurs in our derivation of 
the main result. It is stated, along with the main result, in Section III. 

As in Ref. 1, the complex form of the filter input is 

s(t) = exp [iu t + vp(t)] (1) 

where the carrier frequency is w = 2t/ and the signal is carried by the 
real input phase angle <p(t). Let the filter have the transfer function G(j) 
and the impulse response g(t): 

g(t) = r G(J) exp M) dj, 

G(f) = ( g{t) exp (-tofl dt, to = 2t/, 

where the response g(t) may be complex and may be different from zero 
when t < 0. 

The filter is regarded as a bandpass filter for which G(j) is large only 
near ±/ . Let normalized functions be defined by 

t(j) = ~gUoT' y{t) ~ G(M ' (3) 



(4) 



(5) 
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From the definitions and the Fourier relations (2) it follows that 
y (t) = f" T(/) exp M) df, T(/) = J y(t) exp (-iut) dt, 

I" y(t) dt = 1, r*(-/) = f 7*(0 exp {-i*>t) dt, 

where the asterisk denotes "conjugate complex." 
The filter output corresponding to the input s(t) is 

s o (0 = / g(u)s(t - u) du 

J_oo 

/CO 
g(u) exp [io) (t — u) + i<p(t - u)] du 
-oo 

= G(1o) ] y(u) exp [i<p(t - u)] du\ exp (iu> t) 

= exp (-a - iPo){R(t) exp [id(t)]} exp (ia> t) 

= exp (-a - */3o)!exp [i'0(O]} exp (ico t). 

The definition of y(u) is used in going from the second to the third line. 
In going from the third to the fourth line, the attenuation and phase 
shift, a , /3 at the carrier frequency have been introduced by writing 
G(/o) as exp (-«„ - #o). The complex phase angle 0(Q is related to 
the envelope R(t) exp (-a ) and phase angle -/3 + 0(0 of the output 
[a and /3 are constants which do not depend on <p(t)] by 

exp [i&(t)] = R(t) exp [id(t)], *0(O = In R(t) + i0(«), (6) 

6(t) = Re 0(f), lnB(0 = -Im0(i). 

Comparing the third and fifth lines of equation (5) leads to 

0(0 = -i In f y(u) exp [i<p(t - u)] du\- (7) 

The analysis may be simplified by introducing the "linear portion" 
$(t) of 0(0- Working with the case in which the input <p(t) is small gives 

@(i) = -i in 1 + i J y(u)<p(t - u)du+ • • • J 
~ (— *)t I t(w)<^ - u ) du, 

J-a, 
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and this leads us to define the linear portion of 0(0 as 

*(0 = f y(u)<p(t - u)du. (8) 

J_00 

The complex phase angle 0(0 may be separated into its linear and 
nonlinear portions by adding and subtracting i$(t) in the exponent in 
equation (7): 



0(0 = *(<) - i In 



/ y(u) exp [i(p(t — u) — i$(t)] du • 



(9) 



Adding and subtracting 1 in the integrand and using the fact that the 
integral of y(u) is unity gives the fundamental relation 

0(0 = $(*) - i]n[l +K(t)] (10) 

where 

K(t) = J duy(u){exp [i<p(l - u) - **(<)] - 1). (11) 

In most cases of practical interest, K(t) tends to be small. When 
| K(t) | < 1, expansion of the logarithm gives the series 

0(0 = HD - iK{t) + | K\t) - | K\t) + ■ • • (12) 

upon which our analysis is based. When <p(t) is gaussian, | K(f) | will 
occasionally exceed unity. It appears that results obtained from the 
series of equation (12) represent the first few terms of an asymptotic 
series. This is further discussed in Appendix F. 

If <p(t) is small for all values of t, expansion of the exponential in the 
definition of K(t) [equation (11)] shows that K(t) is 0(<p 2 ). Our main 
expression for the distortion, given in Section III, neglects terms of 
order <p 8 . For this accuracy, equation (12) can be written as 

0(0 = <l>(0 - iK(t) + I K\t) + 0(<A (13) 

Since the variable portion 9(t) of the output phase angle is the real 
part of 0(0, the dc portion, 6 dc , of 0(0 is the average value of Re 0(0- 
When the input <p(t) is a stationary gaussian process with zero mean, 

6 dc = Re (0(O)av 

= Re (-iK{t) + | K 2 (o) av + 0(^ 6 ) (14) 

= Im (K{t) - 2~ 1 K\t))„ + 0(/) 
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where ( )„ denotes "ensemble average" and ($($)» is zero because 
$(0 depends linearly on <p(t). Notice that from equation (5), the total 
output phase angle is 0(0 = -ft, + 0(0 and that the dc part of 0(0 is 

e dc = -&> + e de . 

The two-sided power spectra W$(j), W e (f) o f 0(fli *(0 contain the dc 
spikes (-A, + 0„c) 2 8(J), 6 2 de 8(j), respectively. Furthermore, 

Wi(f) = (# - 2M-«) 8(f) + W,(f). 

Here 5(/) is the unit impulse function. 

In the following work it is convenient to ignore /3 , and we shall 
call 0(0 itself the output phase angle. 

III. STATEMENT OF PRINCIPAL RESULTS 

In all of the results stated here, the input phase angle <p(t) is gaus- 
sian with zero mean. The two-sided power spectrum of tp{t) is W v {f). 
The two-sided power spectrum of the output phase angle 0(t) is W g {f). 

3.1 Second and Third Order Modulation Terms in W e (j) 

The principal result given in this paperf is an expression for Wed) 
which contains all of the second and third order modulation terms: 

W e (j) = 0L 8(f) + WM) I U(ft + U*(-j) | 2 

+ i /I dp w &w& - p) i T (p> / - p) - T *(-p> -/ + p) i 2 

+ ^ /" dp f^ da W v (p)WMW v (j - P - o-) 

• I S(p, a, f - p - a) + S*(-p, -<r, -/ + p + <r) | 2 
+ 0(<p°W,). (15) 



t Note added at press time: Equation (15) gives essentially the first few terms 
of a general expansion due to A. Mircea, Rev. Roum. Sci. Tech. — Electrotechn 
et Energ., 1967, t. 12, No. 3, pp. 359-371, and Proc. IEEE (Correspondence), Octo- 
ber 1966, 54, pp. 1463-1466. I regret the oversight of Mircea's excellent work. Use 
of his results would have substantially improved this article. 
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Here the dc part, d dP , of 0(t) is the imaginary part of 

- \ /_" dp £ da wmwms(-*, -P)[m«, p) - T(<r + p)] 

+ Otf) (16) 

and 

T{p, j - P ) = S( P) j - p) 

+ P da WM[2S(a, p)S(-a, f - p) - S(a, j - a) 

•'—oo 

- T(a)T(-a)S(p, j - p) + S(p + a,f- p- a)] (17) 
U(f) = T(/) + f dp W v (p)V(p)S(-p, f) 

J-to 

+ r dp r daW v (p)W v (a){-±T(p + °)S(-p-a,f) + T(a) 

•1—00 J— 00 

•[3S(-<r, p)S(-p, f) - Sip, j- p-a) + S(p-a,j- p)]| . 

(18) 

The r(/) is the normalized filter transfer function defined by equa- 
tion (3) and the functions S are discussed in Appendix B. They depend 
only on the filter. That is, they are independent of W v (f), and are 
defined symbolically by 

s^ , • • • , O = fl [y xk - rfo)] (19) 

where the power y z of y is to be replaced by T(z) after multiplying 
out the product. The S's are symmetric functions of their arguments. 
For n = 2 and 3, 

S( P , a) = y> + * - y"T(a) - y'T(p) + T(p)T(a) 

= T(p + *)- T(p)T(a), (2Q) 

S(p, a, v) = Tip -f- <r + v) - T(p + a)Yiv) 

- Tip + v)Tia) - Tia + r)T(p) + 2r(p)r(«r)r(v). 
The S(p, a, v) of Ref. 1 is the negative of the one used here. 
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In the "order of" symbol appearing in equation (15) the W v enters 
for dimensional reasons. This is in line with 

f df W,(j) = (At)).. (21) 

J-to 

as can be seen by integrating both sides of equation (15) from / = 
— oo to / = 00. 

In some instances the expression for Wg{j) is useful when rms <p(t) 
is large but rms d<p{t) fdt is small. 

Bedrosian has computed curves showing the second order distor- 
tion for the case of quadratic phase shift, that is, for r(/) of the 
form exp[iA/ 2 ]. 2 

As equation (15) stands, it is oriented towards phase modulation. 
For frequency modulation the time derivatives <p'(t) = d<p(f)/dt, d'(t) = 
dd(t)/dt replace the phase angles <p(t), d(t) as the items of interest. Since 
the power spectrum Wf(j) of 8'(t) is equal to (2irj) 2 W 6 {j) , multiplying 
(15) by (27r/) 2 converts it into an expression for W 8 -(f). On the right side 
of (15), the factor W^H) in the first line is replaced by W,,<(/) and the 
W v 's appearing in the integrands may be transformed into W v -'s with- 
out introducing infinities. The last statement is seen to be true for the 
second order modulation integral when we write 

WMW,(j - P ) = wApWAI - p) WVC/ - p)" 2 

and observe that the product p~\f — p)~ l T(p, j — p) remains finite 
even when p and (/ — p) approach zero. The third order term may be 
treated in a similar way. 

In many applications W v -(j) is proportional to D 2 where 

* - M>„ 

and D is the rms frequency deviation in cycles per second. Then the 
second and third order modulation integrals in (15) are proportional to 
D 4 and D 8 , respectively, as D tends to zero. This suggests that the 
remainder term, 0(cp 6 W v ), is proportional to D 8 . For this reason we shall 
sometimes refer to (15) as the "small deviation" approximation. When 
the FM signal to crosstalk ratio in dB is plotted as a function of log D, 
the behavior of the resulting curve as D — * can be computed from (15). 
Indeed, if the second order modulation predominates, (15) furnishes an 
asymptote to the curve with a slope of 6 dB per octave. If, because of 
symmetry in the filter, the second order modulation term in (15) is 
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zero, the third order term gives an asymptote with a slope of 12 dB 
per octave. 

Some idea of how equation (15) begins to fail as D increases from 
may be obtained by considering the case <p(t) = A sin <a a t. Then the rms 
frequency deviation is D = f a A/2* and the filter output is 

«o(0 = Z Jn(A)G(f + nf a ) exp [i(o> + m> a )t] 

-co 

where J„(A) is a Bessel function and a) a = 2irfa . Consider only the 
second harmonic. It is proportional to J 2 (A), and the approximation 
underlying (15) is roughly equivalent to replacing J 2 (A) by A 2 /8, the 
leading term in its power series. The value of A which makes A 2 /8 
exceed J 2 (A) by 3 dB is A « 2.0 and the corresponding D is 1.4 /„ . If 
the baseband of a gaussian FM wave were flat and extended from to B, 
the expected number of zeros per second would be 1.16 B. This is the 
same as the number of zeros of A sin <a a t with /„ = 0.58 B. This represen- 
tative value of j a leads to the estimate that (15) will be in error by 3 dB 
when D « (1.4) (0.58)B « 0.8 B. Comparison of (15) with experimental 
values indicates that the 3 dB error point typically occurs when D lies 
between B/2 and B. 

3.2 Power Spectrum of ad(t) -{-bin R(t). 

Equation (15) for Wg(f) may be modified to give information re- 
garding the fluctuation of the envelope R(t). This information may 
be of interest, say, in determining the distortion produced by "AM to 
PM conversion." 5 More generally, suppose that one is interested in the 
power spectrum Wo, (/) of 

x(t) = ad(t) + 6 In R(t) = Re [(a + ib)®(t)] (22) 

where a and b are arbitrary real constants. Then W a (f) is given by 
an expression obtained from equation (15) upon replacing U(f), 
T(p,j- P ) and S( P , a, f- P -<r) by (a+ib)U(f), (a+ib)T( P , f- P ), and 
(a+ib)S( P , a, f— P —<r), respectively, so that U* (— /) is replaced by 
(a— ib)JJ* (— /), and so on. (See Appendix E.) 

3.3 Second and Third Order Modulation Terms for "Small and Slow" 
Frequency Deviations 

The expression (15) for W e (f) simplifies when 
(i) T(f) can be expanded as a power series 
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r <fl = ^ = 1 + Sf (23) 

and 

(it) the effective spread of W v (j) is so small that the T's used in 
equations (15) to (20) can be replaced by the first few terms of their 
power series expansion. Roughly, this means that the top baseband 
frequency is small compared with the filter bandwidth. The instantaneous 
frequency changes slowly in comparison with the envelope of the impulse 
response of the filter, and the quasistatic case is approached. With these 
assumptions the resulting simplified form of W e (j) is given by equation 
(126). A more complete form of the small and slow deviation approxi- 
mation is given in equation (133) which brings out the asymptotic 
nature of the results. 

The sum of the second and third order modulation terms given by 
the integrals in equation (126) [which are the simplified versions of 
the corresponding integrals in equation (15) ] is 

W' s (f) = 2- 1 (X 2 ,- + 2~ 1 D% i y r dp W v ( P )W v (f - p)p\j - p) 2 

+ d-'CX.*) 1 f dp r da WMWMW,(j - p - <r)pV(/ - p - v) 2 . 
J-» J-co ( 24) 

Here W l „(j) is the portion of W 9 (j) which gives the interchannel inter- 
ference, that is, the noise a listener would hear in an idle channel in a 
multichannel frequency division multiplex angle-modulation system. 
D is the rms frequency deviation in cps, 



D 2 = 



2tt. 



(25) 



where ( ) av denotes ensemble average. The quantities \ ni are the im- 
aginary parts of the semi-invariants X„ defined by the expansion 

lnr(/)=f;^f- (26) 

71=1 ' t; 

Equations (125) express the first five X„'s in terms of the first five a„'s, 
the coefficients in the expansion of T(j). 

The corresponding approximation for the power spectrum W T x (j) 
of x = a6(t) + b In R(t) [see equation (22)] is obtained by replacing X n , 
in equation (24) by (aX ni + b\„ r ) where X nr denotes the real part of X n . 
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In the important FM case in which the baseband signal <p'(t) has 
a flat power spectrum, the two-sided power spectrum of <p'(t) may be 
taken to be 

WAf) = f°' lU>B (27) 

[(2rD)*/(2B), |/| ^B 

where the baseband extends from / = to / = B (that is, from —B 
to +J5 for two-sided spectra), and D is the rms frequency deviation. 
Then 

w,<fl = WAD/Vvff = j°' ' / ' > B (28) 

[D 2 /(2Bf), \f\*B. 

When this expression is used, the integrals in equation (24) may be 
evaluated and it is found that, for ^ | / | ^ B, 

W',(f) = ^ (2J5 - | / |)(X 2 , + 2- 1 Z> 2 \ 4 ,.) 2 + ^3 (SB 2 - / 2 )(X 3 ,) 2 , 

^ | / | g B. (29) 

The average signal power in an elementary frequency band extending 
from/ to / + A/ in the input base band is W 9 ■ (/) A/ (radians per second) 2 . 
The ratio of the interference power to the output power in the same 

elementary band is 

wUD a/ _ wM m 

w e .(f) a/ w,<j) K6V) 

For the flat baseband FM case we may approximate W g (j) by W v (j) = 
D a /(2Bf) and use equation (29). This leads to the approximation 

Kfr^[( 2 -^ +2 - w+ f( 3 ^M- 

(3D 

for the ratio of the interchannel interference power to the signal power 
in the elementary band (/, / + A/). This ratio has meaning only if 

I / I S £. 
Liou has given an approximation which is equivalent to equation 

(24) for W T g (j) with several more terms included. 6 This approximation 

is discussed in Section XII. 

The small and slow deviation approximation described above gives 
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results which agree well with Monte Carlo computations made by C. L. 
Ruthroff. 6 For illustration we take the simplest of Ruthroff's cases, 
the one in which the transfer admittance of the filter is 

G ® = 1+2, + ^ + ,' - * = * - »"■ • (32) 

Here /„ is the carrier frequency and /„ is the filter semib and width. 
Putting / = /. gives (?(/„) = 1. Putting / = /' + /, gives x = */'//. 
and 

mrCfO-m w + » 



G(f B ) 
= -ln(l + 2x + 2x 2 + x 3 ) 



= _ 2 |L + 2 |__ 48 |_ + 



(33) 



oo //n 

The last line follows from the series of equation (26) defining the Xn's 
as the coefficients in the expansion of In T(J'). In going from line 2 to 
line 3, the logarithm is expanded by setting a Y = 2, a 2 = 4, <x 3 = 6, 
a„ = for n > 3 in 



In (l + Z a-a^/nl) = Z M7w! 



and by using the expressions (125) for the X n 's. 

Substituting x = if JU in line 3 of equation (33) and comparing the 
result with the last fine gives 

X, = -2(*/7 1 ), ^ = 0, X 3 = 2(C) 3 , X 4 = 0. 

Hence X 2 < = 0, X 3 , = —2j~ 3 , X 4 , = 0, and the approximation of equation 
(31) for the ratio of the interchannel interference power to the signal 
power leads to 

- 10 "*■• 31 w - 10 log " Br ( 3 " b)] ■ ^ 

In his Fig. 15 Ruthroff has plotted values of 

-101og 10 [in(/) N/Weij) A/] 
for several different values of D/B and f/B with B = 1 MHz and 



98 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1969 

/„ = 119 MHz. 6 The agreement with our equation (34) is good at 
f/B = 1.0. Our D/B is the same as Ruthroff's a/W and our f/B = 1.0 
corresponds to Ruthroff's slot 10. At f/B = 0.4, equation (34) gives 
values which are about 3 or 4 decibels less than the Monte Carlo values, 
but this may still be regarded as good agreement. 

Similar agreement is found when the small and slow deviation approxi- 
mation is applied to a number of the other cases examined by Ruthroff . 

3.4 A "First Order" Approximation 

Although equation (15) is useful in some FM distortion problems, 
in some cases it is of no help. One example concerns the distortion 
produced by an ideal filter centered on the carrier frequency and having 
a semibandwidth exceeding 3B, where B is the baseband. That is, 
W? (/) is for | / | > B. In this case the distortion is produced by modula- 
tion terms of order higher than the third, and these are neglected in 
equation (15). 

For such problems "first order" approximations can sometimes be 
used. The term "first order" refers to the approximation In (1 + x) £zi x 
where x is of the nature of K(t) in equation (10); it does not refer to the 
order of the modulation products in x. Different choices of x lead to 
different first order approximations. The first order approximation 
given by the first two terms in the series of equation (12) for €>(t) is 

d(t) « Re *(*) + Im K(t). (35) 

The output phase angle 0(0 may be written as the sum 

0(0 = e dc + e ( {t) + e nt (t) 

where 6 { (t) = Re $(0 is the "linear portion" of 0(0, and d nl (t) given by 

MO = 0(0 - t/fi - ^ (36) 

= 0(0 - Re *(/) - 6 de 

is the time varying part of the "nonlinear distortion" in 0(t). The first 
order approximation for 6 nt (t) corresponding to the first order approxi- 
mation of equation (35) for 6(t) is 

d nl (t) « Im K(t) - 6 de « Im K(t) - Im (K(t)) av = y(t) (37) 

where 

y{t) = Im [Kit) - (K(t)U. (38) 

The work of Section VI, which is part of the derivation of equation 
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(15), shows that the power spectrum of y(t) is 

W v (j) = 2" 1 Re [P(j) + Q(j) + P*(-/) + Q*(-f)] (39) 

where 

P(/) = — - / rfT exp ( — tor) / dw / dv 7(11)7 (v) 

•exp [a(u) + a(y)]|exp [2c(u,v, r)] — 1} 

QCf) = o / rfr exp (~* WT ) J rfw J dvy*(u)y(v) 

•exp [a*(w) + a(y)]|exp [2c(u, v, t)] — 1} 

a(K) = -5 f d/ W ¥ (J)H.(j)H u (-f) 

1 J ~'° (40) 

c(w, », r) = -| /_" dj W v (i)H u {-1)H v (f) exp (ur) 

d(u, V, r) = \ J" dj W,(j)H*(j)H.(j) exp (fcr) 

#„(/) = exp (-tort*) - r(/), co = 2t/. 

The function Q(/) is real when / is real, and P(j) is an even function of /. 
The first order approximation W„(J) for the power spectrum of d n( (t) 
contains some higher order modulation terms which are not contained 
in our main equation (15) for W B (f); conversely, equation (15) contains 
terms which are not in W y (J). In using the first order approximation of 
equation (35), which may be rewritten as 

e(t) = Re*(0 + Im K(t) + 0(K 2 ), 

one should guard against throwing away* [in the 0(K 2 ) terms] quantities 
which are of the same order as those being computed from Re <£(0 + 
Im K(t) (the leading term). Although each case requires its own in- 
vestigation, it is helpful to remember that K{t) is 0(<p 2 ). Furthermore, 
when t(0 is real, Im K(t) is 0(<p 3 ). Also when T(j) = 1, K(t) becomes 
0; and when 

| r(/) - 1 | < e « 1 

for all real values of / (as in the case of small wave-guide echoes), it 



*This type of error has been discussed by Enloe, Ruthroff, Gladwin, and 
Medhurst in Refs. 7 and 8. 
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may be conjectured that K(t) itself is 0(e) irrespective of how large 
(<p 2 )„ may be. 

The approximation 6 Ht (t) « y(t), is not quite the same as earher first 
order approximations. 3 • 4,9_12 It is a closer, but more complicated, 
approximation because <p(t - u) — *(<) is used in the integral equation 
(11) for K(t) instead of <p(t — u) — (p(t) as in the earher approximations. 
Appendix D gives some results obtained when the present analysis is 
repeated with <p(t — u) — <p{t) in place of <p(t — u) — $(i). 

Equation (39) for W y (j) gives a first order approximation for the 
power spectrum of d nl (t) = 6(t) - Re *(<) - 6 de . The corresponding 
approximation for the power spectrum of 

In R(t) + lm*(0 - [In R(t)] de tt Re [K(t) - (K(fi)„] = x(t) 

is 

W.(f) = 2" 1 Re [-P0) + Q(j) - P*(-f) + Q*(-/)]. 

3.5 Simplification When Filter has Symmetry T(—f) = T*(f) 
When the filter has the symmetry 

r(-/) - v*(j) (41) 

about the carrier frequency, the even order modulation terms disappear, 
S*(—x 1 , • • • , — x n ) becomes equal to Sfa , •■• , x n ), and equation (15) 
becomes 

W e (j) = W 9 (f) | U(i) | 2 

+ 1 r d P r ^ w^p)wmw,{i - P - *) 

O J-00 •'-00 

■\S(p,cr,f- p-<x) | 2 + 0(^,). (42) 

Here U(J) is still given by equation (18) and S(p, <r, v) by equation (20). 
This expression for W 6 (j) agrees with one of the main results of Ref. 1 
when the double integral in equation (18) for U(j) is assumed to be so 
small that it may be neglected. 

When r(— /) is equal to T*(j), the coefficients a„ in the power series 
of equation (23) for r(/) are real when n is even, and imaginary when 
n is odd. The same is true for the X B 's of equation (26). Hence X 2i , \ i{ 
are zero and the second order modulation terms disappear from the 
small and slow deviation approximations equations (24), (29), and 

(31) for Wl(f). 

The relation r(-/) = T*(f) implies that y(u) is real and that H u (-f) 
is equal to H*(j). Then a(u), c(u, v, t), and c(w, v, t) are real and 6(u,'v, t) 
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is equal to -c(w, v, t). Both P(J) and Q(J) become even and real. Equa- 
tion (39) for the power spectrum of y{t), that is, the first order approxi- 
mation for the power spectrum of the nonlinear distortion 6 nt (t), be- 
comes 

W v (j) = P(/) + Q(f). (43) 

Here the triple integral for P(J) is the same as that given in equation 
(40); and the triple integral for Q(j) may be obtained from the integral 
for P(J) by changing the sign of 2c(w, v, t). Hence equation (43) be- 
comes 

drexp(-^T) / du \ dv y(u)y(v) 

-00 * ""00 " —oo 

•exp [a(u) + a(v)] sinh [2c(u, v, t)] (44) 
where w = 2x/, and a(u), c(u, v, t) are given by equation (40). 

IV. INITIAL EXPRESSION FOR THE POWER SPECTRUM OP 6{t) 

When 0(t) is a stationary noise process its two-sided power spectrum 
W B (J) is the Fourier transform of its autocovariance: 

W B (f) = f exp (-icoT)(d(t)d(t + T)) av dr, u> = 2tt/. (45) 

Denoting functions with arguments t, t + t by subscripts 1, 2 and using 
6(t) = Re 0(<) = 2- 1 [@(t) + 0*(f)] (46) 

gives 

(d(t) 9(t + r)> nv = <M 2 > BV ( 47) 

= 2- 1 Re[(0 1 2 > ttV + (@* 1 2 ) tt v]. 

The procedure of Appendix A and equation (13) for the complex phase 
angle 0(£) lead to 

<010 a )a, = A(T) + A(-T) + 0(^) 
(0*!0 2 )av = B(t) + 5*(-T) + Oi/p) 

where A(t) and B(t) are the ensemble averages 

A( T ) = <$ 1 (2" 1 * 2 - ^ 2 + i2- 1 Kl) - 2~ l K 1 (K 2 - K§)„ 
B( T ) = <*?(2- 1 *, - itf 2 + ^2-7^) + 2- 1 K1(^ 2 - K 2 2 )) av 



(48) 



(49) 
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The remainder terms in equation (48) are 0(cp 6 ) instead of 0(<p 7 ) be- 
cause the ensemble average of an odd order term is zero. 

It also follows from Appendix A that equation (45) for the power 
spectrum of 6(t) goes into 

W e (f) = 2" 1 Re [P(fl + Q(j) + P*(-f) + Q*(-f>] + OtfW ,) (50) 

P(f) = f exp (-iuT)A(r) dr 

(51) 

Q(f) = f exp (-iW)B(r) dr. 

J-oo 

Although the functions P(f), Q(f) used here are not the same as those 
in Section III, they are of the same nature. 

V. CALCULATION OF AVERAGES NEEDED FOR COVARIANCES 

The equations of Section IV show that the value of (Q(t)d(t + r)) av 
depends upon various ensemble averages of products of <!>(£) and K(t). 
When the input phase angle (pit) is gaussian, these averages may be 
computed by using a result proved in Ref. 1. 

Let L be a linear operator which operates on functions of t, and let 

L exp (tut) = exp (mt) lij), a = 2t/. (52) 

Let (p(t) be a stationary gaussian process with two-sided power spectrum 
W 9 (j). Then 



(exp [iL(p(t)]) av = exp 



-| /_" df w v mm-i) 



(53) 



Setting xL(p(t) for Lcp(t) and comparing coefficients of x 2 in the power 
series expansions of the two sides of equation (53) shows that 

fji w,wm-fi = ([L(p(t)f) av . 

That (exp [iL(p(t)]) av is equal to exp {— 2 _1 ([L^(0] 2 ) Q v} follows from 
the fact that the real and imaginary parts of Ltpify are correlated gaussian 
processes. 

In dealing with K(t) it is convenient to introduce the function J(v, t) 
defined by 

J(v, t) - exp [vp(t + t - v) - &(t + r)]. (54) 

The dependence of J(v, t) on t is ignored because the right side of 
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equation (54) is a stationary random process, and J(v, t) will be used 
only to calculate ensemble averages. The examples, which follow from 
the definition of equation (11) of K(t), 

</Q, lv = (K(t)) av = j\luy(u)(J(u, o) - 1>„ 

{KJU\ V = (K(t)K(t + r)>„, (g5) 

= f du r dvy(u)y(v)([J(u, o) - l][J(v, r) - 1])„ 

J —CO "—00 

(K* X K,)„ = f^ du r dvy*(u)y(v)([J*(u, o) - l][J(v, r) - 1]>„ 

■/— eo "—00 

show that averages of the type (J(u, o)) av , (J(u, o)J(v, r)) av , and 
(J*(u, o)J(v, r)) av are needed. 

To calculate {J(u, o)) av from the general result, equation (53), let 
L be the operator which carries <p(t) into <p(t — u) — $(<)• Replacing 
<£>(£) by the integral which defines it gives 

/oo 
ds y(s)(f(t — s) . 
-00 

The function /(/) associated with L is obtained by setting exp (uat) 
in place of <p(t)\ 

exp (iut)t(f) = /.[exp (tut)] 

= exp [i(x)(i — u)] — I ds y(s) exp [iw(l — s)], 

Kj) = expi-iou) - T(f) -J5T.0). 
Then equation (53) gives 

(exp [*I^(0]>»v = (exp [i<p(t - u) - *'*(<)]>., 

= (/(«, o)) nv = exp [a(u)\ (56) 

where 

a(u) = -~ r dj W v (f)Hu(Mu(-1). (57) 

& J-00 

The functions a(u) and 

ff.fl) = exp (-torfu) - T(j) (58) 



(59) 
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play important parts in the analysis. The present H u (f) is the negative 
of the one used in Reference 1, a change made to simplify the analysis. 
The calculation of (J(u, o)J(v, t)) &v proceeds in much the same way. 
Let 

L<p{i) = <p(t - u) - $(«) +<p(t+ t -v) - $(< + t), 

t(fl = #„(/) + exp (zW)#,(/), 

= -\ /_" df W.(MH„(f)Hui-f) + H.(f)H.(~f) 

+ 2 exp (icor)H u (-j)HM 
= a(u) -f- a(v) + 2c(u, v, t) 
where W v (j) is an even function of / and c(u, v, t) is the integral 

c(u, v, t) = -| f df W v (f)H«(-f)H v (j) exp (i2ir/r). (60) 

Consequently, 

(J(u, o)J(v, T)) av = exp [a(u) + a(v) + 2c(u, v, t)]. (61) 

Similarly, to calculate (J*(u, o)J(v, r)) ttV let 

Lcp(t) = -<p(t -U) + $*(t) + <p(t + T - V) - $(t + r), (62) 

£(/) = —exp ( — zW) + / dsy*(s) exp ( — iois) + exp (tW)i/„(/) 

= -#*(-/) + exp (*W)tf,(/) 

where the Fourier transform of y*(s) is T*(— /). The work of equations 
(59) and (60) goes through much as before with —#*(—/) in place of 
HJj). The result is 

(J*(u, o)J(v, r)> av = exp [a*(u) + a(v) + 2c(u, v, t)] (63) 

where 

c(u, v, r) - | /* d/ W v (f)H*«(f)H v (1) exp (*2tt/t). (64) 

All of the averages needed are given in Tables I and II. Items 1, 3, 
and 6 in Table I have just been computed, and the others may be ob- 
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tained in a similar manner. The entries in the last column of Table I 
may be verified by expressing the a's and c's as ensemble averages 
[see equation (67)] and using (exp [iL<p(t)]) av = exp {— 2~ l {[Up{t)f) av \. 
Table II gives averages of products in which one factor is <p(t — u). 
The first average, (<p(t — u)<p(t + r — v)) av , is the Fourier transform of 
W P (J). The second average, {<p(t — u)J(v, r)) av , is the coefficient of ix 
in the expansion of exp ([iL<p(t)]) where 

Lip{t) = x<p{t - u) + #>(£ + r - v) - *(* + r) _ 

/(/) = x exp (— icau) + exp (icoT)H v (f). 

The third average may be computed in a similar way. 

The following list brings together the integrals a{v), c(u, v, t), 
which appear in Tables I and II: 

a(u) = -| /" dj W v (j)H u {f)H u {-f) 

c(u, v, r) = -\ ^ df W v (f)H«(-f)H*(.f) exp (i2x/r) 

& J-* (66) 

a*(u) = -| J" df Wr(f)H*(j)H*(-I) 

c(u, v, r) = \ |" dj W^j)H*»(f)H,(f) exp (*2x/r) 

where H u {f) = exp (-i2irfu) - T(J), and replacing H u (j) by -H* u {-j) 
in a(u), c(u, v, t) gives a*(u), c(u, v, t). Also 

a(u) = c(u, u, o), c(w, v, 6) = c(v, w, o), 

c(u, v, -t) = c(v, u, t), 6(u, v, -t) = c*(v, u, t), , q „ 

c(u, v, r) = -h([<p(t ~ u) - *{t)I<p(t + t -»)-*(* + r)]) av 

c(u, v, r) = \{[<p(t - u) - $*(t)][<p(t + r - v) - $(i + r)]) av . 

When r(— /) is equal to r*(/), both 7(0 and <£(£) are real; and it 
follows that a(u), c(u, v, t), c(u, v, t) are also real. Furthermore, H u (—f) 
= H*(f) and c(u, v, t) = — c(u, v, t). 

VI. THE POWER SPECTRUM OF K(t) 

The dc portion of the complex random process K(t) defined by the 
integral in equation (11) is the complex constant 
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<K(0>.v = <^,) av = f M du y(u)(J(u, o) - l> av 

J -« (68) 

/CO 
du y(u) { exp [a(u)] — 1 } . 
-00 

This follows from equation (55) for (K 1 ) av and the expression for 
(J(u, o)) av given in Table I. 

The power spectrum of K(t) is the Fourier transform of (K* X K 2 )„ • 
The integral for (K*K 2 ) av given by equation (55) and the ensemble 
averages of J*(u, o), J(v, t), and J*(u, o)J{v, t) given in Table I lead to 

/CO /*co 

du \ dv y*(u)y{v) exp [a*(u) + a(v)] 
-CO " — CO 

•{exp[2c(w,y, t)] - 1}. (69) 

Integrals of the type appearing in equations (68) and (69) may be 
expressed as infinite series involving the S functions [which depend 
only on r(/)] described in Appendix B and the more complicated func- 
tions S n described in Appendix C. Only the first few terms need be 
considered when most of the distortion arises from second and third 
order modulation. 

The definition [equation (167)] of the complex constant S and its 
series expansion [equation (171)] give 

(K(t)) av = S - 1 

= -| /" dp W v (p)S(p, -p) 

+ 1 f° dp r da wMWMSb, *, - P , -*) + o(A 

O J -co J -co CJQ\ 

Expanding exp [2c(u, v, t)] in equation (69) in powers of 2c(u, v, t), 
replacing each c(u, v, t) by its defining integral [equation (66)] with p in 
place of /, and integrating with respect to u and v with the help of 



S-(Pi , • • • > P») = P dv 1® ex P [a(y) ^ II H.( Pk ) (71) 

J-to * = 1 



leads to 



<K*K 2 ) av = i (k^ \*+thr d Pl --- r dpiti w 9 ( Pk )\ 

•exp [i2TT( Pl + • • • + p„)]5*( Pj , • • • , Pn )S n ( Pl , • • • , p n ). (72) 
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When this expression for <i^^K 2 ) ov is put in the integral 

W K (f) = r exp (-iurXKW^ dr, eo = 2tt/ (73) 

J-oo 

for the power spectrum, W K (j), of K(t) and the Fourier transform of 
unity denoted by S(j), the result is 

W K (f) = | (K(t)) av | 2 5(f) + P dr exp (-icor) 

J_oo 

/OO /«0O 

dw / dv y*(u)y(v) exp [a*(w)+a(i>)] { exp [26(u, v, t)] — 1 } 

= | So - 1 | 2 W)+ £i f d Pl ••• /""dp,, 

n _l /tl J -oo •'—00 

■ «(/ - P, Pn)[l[ W»(A)j I SniPl , ' • • , P„) | 2 

= | £ - 1 | 2 8(f) + *%(/) | £,(/) | 2 
+ \ /_" <*P W 9 tp)W,(f - p) | S 2 (p, / - p) | 2 

+ | £ dp f_ M da W v (p)W v (ajWM - p-a) 

■ I S 3 (p, cr, f - p - a) |" + • • ■ . (74) 

The leading terms in the series for S Q , S t , S 2 , £ 8 m terms of un- 
subscripted S's are given by equations at the end of Appendix C. The 
inequality for S„ given in Appendix C may be used to show that the last 
series in equation (74) converges when W v (f) remains finite for all 
values of / and ([<p(t)] 2 ) av is finite. The convergence of similar series which 
will be encountered later will be tacitly assumed. 

VII. "FIRST ORDER" APPROXIMATION FOR POWER SPECTRUM OF 0(t) 

Before taking up the problem of computing W$(f) from d(t) = Re ©(t) 
and 

®(tj = *(*) - iK(t) + I K\t) + 0(<p 6 ), 

which is the same as equation (13), we shall go through a similar, but 
simpler, calculation using the "first order" approximation 

0.00 = 4>(0 - iK(t) +G(*> 4 ). - .-.;..■ •-...■ (75) 
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Neglecting 0(K 2 ) terms in equation (14) for de , and in the covariance 
expressions given in Section IV, shows that 

6 de = Im (K(O)av + 0(* 4 ) (76) 

W B (f) = 2" 1 Re [P(1) + Q(f) + P*(-f) + Q*(-f)] + OtfW v ) 

where P(f), Q(f) are the respective Fourier transforms of 

A(r) = <$ 1 (2" 1 $ 2 - iK 2 ) - 2~ 1 K 1 K 2 )„ (?7) 

B( T ) = (^(2 _I $ 2 - iK 2 ) - 2~ 1 K*K 2 ) av 

Expressions for (K^-K^)*, have been obtained in the preceding sec- 
tion. Repeating the work with K x in place of K* brings in y(u), a(u), 
c(u, v, t), -H u (-p), and (-) n £„(-pi , • • • , -p„) in place of y*(u), 
a*(u), 6(u, v, t), and H*(p), £*(p, , • • • , p„). The result is 

(KJU)^ = [(K^f + r du f X dvy(u)y(v) exp [a(w) + a(v)] 

J— CO "—CO 

•fexp[2c(w,v, t)] - 1} 

= [<* i>.j 2 + Z ^r dpi • • ■ / <H II w,<a) 

„_! 71! J-oo •'-co Lft-1 J 

• exp [i'2xt(p! + ■ • • + Pn)] 

•SnC-Pi , ••• , -Pn)S n ( Pl , ■■■ , P n). (78) 

The remaining portion of A (t) in equation (77) is 
<$ 1 (2- 1 $ 2 - iK 2 )) av 

= du I dv y(u)y(v)(2~ 1 (p(t - u)<p{t + t - v) — i<p(t — u)J(v, t))„ v 

•J— CO "—CO 

= / du J dv y(u)y(v) 

J -to •'-co 

• r d] W r (j) exp [uc(t + w)]|2- J exp [-itov] + #„(/) exp [a(w)]} 

J -CO ( 79 ) 

where co = 2irf and Table II has been used in going from the first equa- 
tion to the second. Integration with respect to u brings in T(— /), and 
integration with respect to v brings in both r(/) and the function <S t (/) 
of Appendix C. 

(* 1 (2' 1 $ 2 - iK 2 )) av = [" dj W v (j) exp (icor)r(-/)[2- 1 r(/) + SM- 
J - a (80) 

Replacing $, , y(u) by $f > y*( u ) causes T*(J) to appear in place of 
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r(— /) and shows that the remaining portion of B(r) in equation (77) 



IS 



<$t(2- l * 2 - iK 2 )) av = r df W,(j) exp (fcor)r*(/)[2- x r(/) + SM- 

The function A(j) is the sum of — 2~ 1 (K i K 2 ), obtained from equation 
(78), and (80). Its Fourier transform is 

P(f) = T^,(/)r(-/)[2- 1 r(/) + sm - 2- 1 [<-^ 1 >„v] 2 s(j) 

— - / dr exp (— i2irjr) j du j dvy(u)y(v) 

~ «/— oo **— OO •'—00 

•exp [a(u) + a(t>)]{exp [2c(u, v, r)] - 1} (82) 

where the leading term follows immediately from equation (80) and 
the Fourier integral theorem. 

The function B(t) is the sum of 2~ 1 (KfK 2 ) av , obtained from equa- 
tion (69), and (81). Its Fourier transform is 

Q(f) = WvW^m-'ni) + sm + 2- 1 1 </q 11v 1 2 «o> 

+ - / dr exp ( — i2tJt) j du j dvy*(u)y(v) 

• exp [a*(u) + a(v)} { exp [2c(w, v, t)] - 1} . (83) 

A first order approximation for W e (J) may be obtained by combining 
equations (76), (82), and (83). Deleting the terms multiplied by 8(f) 
and W v (j) gives the first order approximation to the power spectrum of 
the nonlinear distortion n , (t). This approximation is stated by equations 
(38), (39), and (40) in the section describing the results. We now proceed 
to express the first order approximation for W g (j) as the series given by 
equation (90). 

When equation (82) for P(f) is added to equation (83) for Q(f) and 
the triple integrals replaced by their series, namely, the Fourier trans- 
forms of the series appearing in equations (78) and (72), the result is 

P(f) + Q(i) = w v (f)[V(-f) + r*(/)] 

■ p-Tfl) + sm + 2" 1 (| <*!>„ | 2 - [<*!>„]■) *(/) 



+ \ th r dp, ■ • . r d Pn 8d- P1 — Pn )\fi w. 



U) 



(84) 
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Changing the signs of / and the variables of integration p t , • • ■ , p„ , 
and then taking the conjugate complex shows that the series in the 
expression for P*(— /) + Q*(— /) differs from equation (84) only in 
that the S„ factors are replaced by 

S*(- Pl , ■■■ , -Pn)[-(-) n S*( Pi , ■ ■ ■ , p„) + S n (- Pl , ■•■ , - Pn )]. 

(85) 

Taking (— ) n_1 out of the square brackets in equation (85) and then 
adding the series term in P(/) + Q(f) to the series term in P*(— /) + 
Q*(-/) gives 

I Z n f rfp. • • • f d Pn 8(1 - p, p n )[ll W v ( Pk ) 1 

Z „=i Til J-oo J-oo La.- = 1 J 

• I S.( Pl , • • • , p.) - (-)•«(-/», , • • • , -Pn) I 2 (86) 

for the series term in P(f) + Q(j) + P*(-/) + Q*(-/). 
The term for n = 1 in the series of equation (86) is 

**%</) I S X (J) + #(-/) | 2 . (87) 

From the first line in equation (84), the sum of the other terms in 
P(f) + Q(f) + P*(~f) + Q*(~f) containing the factor W v (f) is 

w r (mn-t) + r*(/>] 

■ p-'rfl) + 2-'r*(-/) + s 1 (j) + £*(-/)]. (88) 

The real part of the sum of equations (87) and (88) may be written as 

WM) I r(/) + r*(-/) + s x (j) + s*(-j)\ 2 (89) 

These results and equation (76) for W e (f) lead to 

w,(j) = el 8(f) + 4-tf,(/) | r(/) + s l (fl + r*(-/) + s?(-fl | 2 

+ 4" 1 f, h f d Pl • • ■ f dp„ 5(/ - Pl Pn )[ fl Tf,(p,)l 

• I S„( Pl , • • • , p„) - (-)"«(- Pi , • • • , -Pn) I 2 + OtfWJ. (90) 

The remainder in equation (90) for W e (j) is 0(<p*W v ) while the one 
in the main result, that is, equation (15), is 0(<p°W 9 ). The result of 
neglecting all 0(ip*W v ) terms in equation (90) agrees with the result 
obtained by neglecting the 0(<p 4 W v ) terms in the main result. This may 
be verified with the help of 
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s 1 a) = -\ r dp w v (p)s(p, -p, f) + ok*) 

2 J -« (91) 

S 2 (p, v) = S( P , v) + 0(<p 2 ) 

which follow from the equations at the end of Appendix C. 

viii. "second order" approximation for the power spectrum of 6(t) 

In Section VII the "first order" approximation to Ws(f) is com- 
puted using the approximation 

®(t) = *(t) - iK(t) + 0(<* 4 ) (92) 

for the complex phase angle %{t). In this section the "second order" 
approximation to W g (f) will be computed using the approximation 
given by equation (13), 

0(0 = $(tj - iK(t) + \ K\t) + Ofe> 6 ). (93) 

The equations needed are given in Section IV. Portions of the ensemble 
averages A (t), B(t) defined by equation (49) have already been obtained 
in Sections VI and VII. The remaining portions needed are 

(fi-'MX , {2~ l KJ^„ , (94) 

for A(t) and 

{i2-'$* x Kl) av , (-2- 1 K* 1 K 2 2 ) av (95) 

for B(t). 

From Table II, 

W^KD^ = «2 _1 f du\ dv f dw y(u)y(v)y(w) 

J— oo J — oo « — 00 

■(<p(t - u)[J(v, t)J(w, t) - J(v, r) - J(w, t) + 1]>„ 

/CO *»00 /»CO 

du I dv I dw 7(^)7(^)7(10) 

'/I df W *® exp Mt + uKVHMlVWz - V], 

co = 2tt/. (96) 

In going from the first to the second equation, symmetry in v and w 
has been used to replace H v (/) + H w (/) by 2H, (/) and we have introduced 
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part of the notation 

U = exp [o(«)] f V = exp [a(v)] } W = exp [a(w)] 

x = exp [2c(u, v, t)], t/ = exp [2c(u, w, t)], z = exp [2c(w, v, o)] 

A = exp [2c(u, v, t)], ij = exp [2c(u, w, t)], U* = exp [a*(w)]. 

(97) 

As in equation (80), integration with respect to u brings in r(— /), 
and integration with respect to v and w brings in the functions &i (/;)> 
£,(/) of Appendix C: 

(j,2-^ x K\) av = -T df exp (io>r)W ,(/)r(-/)[S 10 (/;) - £,(/)]. (98) 

J_oo 

The corresponding portion of B(t), ^'^Kl)^ , is equal to the 
expressions obtained when y(u) and T(—f) are replaced by y*(u) 
and r*(/) in the right sides of equations (96) and (98). 

The last portion of A (t) is 

/CO /iCO i»QO 

du dv dw y(u)y(v)y(w) 

•<[J(W, o) - 1][J(«, r) - 1][J(«;, r) - 1])„ 

= C t + A(r) (99) 

where Cx is independent of r and represents the value of equation (99) 
at t = oo . With the help of Table I and the notation denned in equation 
(97), the ensemble average in the integrand may be written as 

UVWxyz - UVx - UWy - VWz + U + V + W-1. (100) 

The only variables in this expression which contain t are x and y. When 
T — > oo , c(u, v, t) tends to and x and y'tend to 1. Therefore the portion 
of equation (100) which contributes to C\ is 

UVWz - UV - UW - VWz +C/ + 7 + TT-1 

and the portion contributing to D x (t) is the remainder 

UVWz(xy - 1) - UV(x - 1) - UW(y - 1). (101) 

The portion contributing to C\ will be ignored since the Fourier trans- 
form of C\ , namely C t 8(f), is part of d 2 dc 8(j), and d de will be treated by 
itself. 
When xy — 1 is written as (x — l)(y — 1) + (x — 1) + (y — 1) 
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equation (101) becomes 

UVWz(x - l)(y - 1) + UV(Wz - l)(x - 1) + UW(Vz - \){y - 1). 

The symmetry in v and w allows the last summand to be replaced by 
the second and hence 



A(t) = / (hi / dv J dw y(u)y(v)y(w) 



•[2~ l UVWz(x - l)(y - 1) + UV(Wz - l)(x - 1)]. (102) 

Expanding (x — 1), (y — 1) in powers of c(u, v, t), c(u, iv, t), respec- 
tively, and integrating termwise in much the same way as in the passage 
from equation (69) to equation (72) leads to 



D 1 (r) = ±^fjy Pl ...f_ 



dp n 



II w. 



I Pa- J 



• exp [i2 7 rr(p 1 + • • • + p„)] 

•S n (— pj , • • • , — p n )[£„o(pi , • ■ • , p„;) — Snip,. , • • • , p»)] 

■/I dai " ' .C rf<r '" [ft Tl ^ (pi) ][ft tf >*)] 

• exp [i2iTT(p s + • • • + p„ + a, + • • • + O] 

•«Sfn+m(— Pl , • • • , — Pn , —0"l , ' • • , —Cm) 

■sUpi , ••• ,Pn ;^i , ••• ,0- (103) 

When the last portion of B(r) is written as 

(-2- l K*K 2 2 ) av = C 2 + D 2 (r) (104) 

the work goes through much as for C x + D x (t). The functions y*(u), 
J*(u, o), U*, £, and y replace y(u), J(u, o), U, x, and y, respectively. 
The functions H u ( — p h ), and H u { — <j t ) in c(u, v, t), c(u, w, t) are replaced 
by —H*(p k ), and —H*{(T t ). This carries x, y into £, y and causes £„(pi , 

• • • , — Pn) to be replaced by (— ) n £*(pi , • ■ ■ , p„). A similar replacement 
holds for S n+m . 

The resulting expression for D 2 (t) is obtained by changing the sign 
(because —K.% replaces K,) of the expression (103) for Di(t), and then 
replacing S n {-p x , • • • -p„) and >S n+m (-p l , • • • , -<r m ) by (— )"/St(pi , 

• • • , p n ) and ( — ) n+ '"*S>* +m (pi , • • ■ , <r m ), respectively. 
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Now that expressions for the portions (94) and (95) of A(t) and B{t) 
have been obtained (in effect), there remain two problems: 

(i) taking their Fourier transforms to get their contributions to P(/) 
and Q(f), and 

(«) adding these to the first order approximation for W t (f) given by 
equation (90). 

The Fourier transform of (pT^^KD^ follows from equation (98) 
and the Fourier integral theorem. The Fourier transform of (i2 _1 <f>*.K;;)„ v 
may be obtained in much the same way and consequently the contribu- 
tion of these terms to P(j) + Q(j) is 

(" dr exp (-te»T)[<i2~ l *i.K2>., + (»2" 1 ** 1 Kj>.v] 

•'-00 

= -W v (f)[T(-f) + r*(/)][<S,„(/;) - SM- (105) 
Consequently their contribution to the right side of 

We(f) = 2" 1 Re [P(D + Q(1) + P*(-/) + Q*(~m + 0( v *W r ) (106) 
[from equation (50)] is 

2~ 1 W,(D Re [r(-/, 4- r*(/)][S,(/) - s I0 (/;) + si(-/) - Sf (-/;)]. 

(107) 

The Fourier transform of the portion Dj(t) of (2~ 1 K 1 Kl) av is obtained 
by replacing exp [t2irr(pi + ■ ■ • + P„)] and exp [i2itt(pi + ■ ■ ■ + O] 
in equation (103) for D^t) by «(/ - p, - ■ ■ • - P») and 5(/ - p, - 
... _ - m ) | respectively. The Fourier transform of D 2 (t), from 
(— 2~ 1 K^K2) av , can be obtained similarly. The sum of these two Fourier 
transforms gives the contribution of AM 4- D 2 (t) to P(f) + Q({)- 
Changing the signs of / and the variables of integration Pi , • • • , <r m , 
and then taking the conjugate complex, gives the contribution of 
D^ T ) 4. £) 2 ( r ) to P*(-j) + Q*(-/). When the two contributions are 
added, it is found that the contribution of DM 4- D 2 (r) to P(J) + 
Q(/) + P*(-/) 4- Q*(~i) is 



E -i / d Px ■ • ■ / dp. n ^(pO 

■*(/- Pi - ••• " p-)[S*-(pi , • 
■[S n ( Pl ,•••)- (-) b «(-p, , • 

- SMi , •■•) + (-)"^*n(-p 1 , 



■) - (-) n 5 n (-p, ,•••)] 



■) 



■)] 
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1 °° _°o I poo /»O0 

[n ^(^)][n w f M\ «(/ - pi o 

■wa..(A, •••)-(-r m ^ +m (-Pi, •••)] 

■LSUGh ,■•■)- (-) n+m SL(- Pl ,-■■)] (108) 

where the complete arguments of the S functions are shown in equation 
(103). 

The desired expression corresponding to equation (106) for W g (j) 
is obtained by adding the significant terms in the first order approxi- 
mation of equation (90) for W g (f) to the second order terms given by 
equations (107) and (108). The remainder term, 0(<p 4 W v ), in equation 
(90) can be ignored because the significant terms are obtained from 
$(£) — iK(t) without approximation [compare equations (92) and (93) 
for 0(0]. The result is 

w e (f) = e 2 dc 8(f) + 4r l w 9 (f) | r(/) + s l (f) + r*(-/) + #*(-/) | 2 

+ 4" 1 Z^ f dp, ■■■ P d Pn 8(f- Pl Pn)[n^(P*)l 

n=2 ni J-„ J-„ Lt-1 J 

■|s»(pi, ••• ,p n ) - (-ys* n (- Pl ■■■ - Pn ) | 2 

+ expression (107) + 2" 1 Re[expression (108)] + 0(<p°W v ). (109) 

The next section is concerned with the elimination of all 0(<p*W v ) 

terms from the significant portion of equation (109). When these terms 
are eliminated, the result is the "main result" stated in equation (15). 

IX. ELIMINATION OF HIGHER ORDER MODULATION TERMS FROM W${f) 

In this section all terms of 0(<p a W v ) in equation (109) for W e (j) 
will be discarded, that is modulation terms of order higher than three 
will be discarded. Since the integral of W v (f) is 0(<p 2 ), all terms in 
equation (109) containing the product of four or more W„'s may be 
dropped immediately. 

First consider the terms which explicitly contain the product of three 
W<p's. This corresponds to n = 3 in the single series in expressions (108) 
and (109), and to the pairs of values n = 1, m = 2; n — 2, n — 1 in 
the double series. The contribution of the double series can be dis- 
carded because it is 0(f / W\<p 2 ) = 0(<p 6 W v ), the functions S 21 and 
S ia being 0(<p 2 ) [from S(f) = 0, £,(/) = 0(<p 2 ) and Appendix C]. The 
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n = 3 term in expression (108) can also be discarded because, from 
Appendix C, S 3 — S 30 is 0(<p 2 ). Using S s — S = 0(<p 2 ) in the n = 3 
term in equation (109) shows that the contribution to W e (J) of the 
terms which explicitly contain the product of three W v 's is 



24 il dp /_" d(r W M W M W *V) 



• | S(p, <r, v) + S*(-p, -<r, -v) | 2 + OtfW,) (110) 

where v = / — p — <r. 

Next consider terms which explicitly contain the product of two 
Wy's, namely the terms n = 2 in the single series and n = m = 1 in 
the double series. When we put p x = p, p 2 = / — p = v in the single 
series terms, and p t = p, <x x = f — p = v in the double series term, all 
of the integrands contain the factor 

P* = [S%( P , v) - S 2 (- P , -v)] 
and the contribution of their sum to W g (j) can be written as 

I /." dp W *bW,<?) Re [fi*(fi + 2 7 )] 

where /3 is 0(1), and y is not the earlier y(u). Here 
7 = S a (p, v) - S 20 (p, v;) - S ll (p; v) - S*(- P> -v) 

+ S* (- P , -v-) + S*(- P ; -u) 
is 0(<p 2 ) since both S 2 — S 20 and S u are 0(<p 2 ). Furthermore, 

Re [0*03 + 2 7 )] = | /3 + 7 I 2 - 77* = I ^ + 7 I 2 + 0(A 

/3 + 7 = t{p,v) - f*(-p, -v), (111) 

f(p, v) = 2S 2 (p, „) - S 20 (p, v]) - S n (p, v). 

The equations at the end of Appendix C may be used to show that 
f(p, v) is equal to T(p, v) + 0(tp) where 

T(p, v) = S(p, v) + f da WMl-hS(<r, -a, p, v) 

J-eo 

+ m*> -*)S(p, v) + S(a, p)S(-a, ,)] (112) 

and consequently the contribution to W $ (j) of the terms which explicitly 
contain the product of two WpB.is '. - . - - 
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l /_! dp W M W *® I r <* ') " T *(-P> -") I' + OtfWJ (113) 
where v — f — p. 

Now consider the terms in W B (j) which are multiplied by W V (J). From 
equations (109), (107), and the term for n = 1 in the single series in 
equation (108), the sum of these terms is W v (j) multiplied by 

4" 1 1 m + &(/) + r*(-/) + s*(-j) | 2 
+ 2" 1 Re [r(-/) + r*(/)][S,(/) - 5 10 (/;) + S*(-f) - S&(-/0] 
+ 2- 1 Re [Si® + £,(-/)][£,(/) - 5 10 (fO + S1(-fl - Sfc(-/;)] 

= 4" 1 | a + j8 | 2 + 2" 1 Re (a* + /3*)( 7 ) 

= 4~ 1 [(a + 0)(a* + jS*) + (a* + /3*) 7 + (a + /3) 7 * + 77* - 77*] 

= 4" 1 | a + + 7 T ~ 4 _, 77* 

= 4- | £/(/) + #*(-/) I 2 - 4 _ Vy*. (114) 

Here, with /3 and 7 different from those in equation (111), 

a = r(f) + r*(-/) = 0(1) 

= s 1 (f) + SK-f) = o(<p 2 ) (U5) 

7 = Sx(/) - &„(/;) + stf-fl - SU-f-) = Otf) 
t)(1) = r(/) + 2£ 1 (/) - £ 10 (/;). 

The equations at the end of Appendix C may be used to show that 
#(/) is equal to U(f) + 0(<p°) where 

U(f) = T(f) -\f dp WMS(p, -p, f) + P dp r da W v (p)WM 

— J— oo •'—00 «—oo 

•[|£(p, <r, -p, -o-, /) - jS(p, —p, f)S(a, —a) 

- ±S(p, f)S(a, -a, -p) - |S(p, a, /)£(-p, -a-)]. (116) 

Since 77* is 0(<p 8 ), the terms in W e (j) which are multiplied by W v {j) 
can be written as 

W,(j) \U(j) 4- t/*(-/)| 2 + 0(/lF,). (117) 

Finally consider the dc spike, 6 2 de 8(j), in W B (j). From equation (14), 
the dc component of 0(0 is 

d dc = Im (Kit) - 2~ 1 K 2 (f)) liV + 0(<p 6 ). (118) 
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The value of (K(t)) av is given by equation (70) and from equation (78) 
with t = 0, 

(K 2 (t)) av = (K(t))l - r dp W v (p)S 1 (-p)S 1 (p) 

J -to 

+ 5 f dp r da W,(P)WMS 2 (-P, -o)S 2 (p, a) + 0(<A 

Z J -co J — oo 

Since Si(p) is 0(<p 2 ), the single integral is 0(<p 6 ) and may be included in 
the remainder term. Squaring the leading term in equation (70) to get 
(K(t)) 2 , combining terms, and using S 2 (p, <j) = S(p, a) + 0(<p ) 
leads to 

(K(t) - 2-7v 2 (0>„v = D e 

D = -\ f dp W,(p)S( P , -p) + f^ dp |_" do WJJiWjp) 

•[*S(p, *, -p. -«0 - *5( P> -p)S(<r, -<r) - KS( P , <r)£(-p, -cr)] 

+ 0(«A (119) 

The imaginary part of D e gives 6 de . 

Addition of equations (110), (113), and (117) shows that 

We(i) = dl 8(1) + W v (f) | U(f) + U*(-f) | 2 

+ i F dp w m w & - p) i T(p ' 1 - p) " r * (_p ' _/ + p)|2 

+ ^ /^ ^ /" d* w v ( P )W v (*WM 

■ | S(p, a, u) + S*(-P, -*, -») | 2 + 0(/TF,) (120) 

where v = f — p — a. This is the same as equation (15) in the state- 
ment of results. However, the expressions for D c , T( P , f — p), and 11(1) 
given in Section III are simpler than the ones given in this section. The 
method of obtaining the simpler expressions will be outlined in Sec- 
tion X. 

X. SIMPLIFIED EXPRESSIONS FOR 6 dc , U (f) , AND T(p, v) 

The expressions obtained for Q dC) U(f), and Tip, v) in Section IX 
may be put in forms better suited to calculation by writing the higher 
order S functions in terms of S functions of two arguments, 

S(p,a) = T(p + <r) - r(p)r(a). (121) 
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These simplified forms are the ones stated in equations (16), (17), and 
(18). 

Since no really satisfactory procedure of reduction was found, the 
expressions given here may not be the simplest. The procedure is il- 
lustrated for the double integral 



= r dp r da wmwm 

J —cc J— en 



•[S(p, -p, cr, -cr) - S(p, -p)S(a, -cr) - 2S(p, cr)S(-p, -<r)] 

which appears in equation (119) for de = Im D c . 

After some cancellation, the general equation (159) for S(p, cr, v, p.) 
shown in Appendix B gives 

Sip, ~P, *, -cr) = 1- (p)(-p) - (v)(-cr) + (p + cr)(-p)(-cr) 

+ (P - o)(-p)(«r) + (-P + «r)(p)(-«r) 

+ (~P - <0(p)(<r) - 3(p)(-p)(<r)(-<r). 

Here T(x) has been written as (x) and (0) = r(0) = 1 has been used. 
When this expression is multiplied by W„,(p)W ip (<r) and integrated 
with respect to p and cr, changes in the variables of integration show that 
the value of / is unchanged by the substitution 

S(p, -p, cr, -o") -» 1 - 2(p)(-p) + 4(p + cr)(-p)(-<r) 

~ 3(p)(-p)(cr)(-«r). 

Here the arrow means "may be replaced in the double integral by". 
Similarly, 

-S{ P) -p)Sia, -<r) - -[1 - (p)(-p)][l - (o)(-*)] 

-> -1 + 2(p)(-p) - (p)(-p)W(-<r) 

-2S(p, <r)S(-p, -a) = -2[(p + cr) - (p)(o)][(- p - a) - (-p)(-<r)] 

-> -2(p + a)(-p - <r) + 4(p + cr)(-p)(-<r) - 2(p)(-p)(<0(-<r). 

Addition shows that the quantity within the square brackets in the 
integrand of I may be replaced by 

8(p + a)i-p)(-<r) - 2(p + cr)(-p - cr) - 6(p)(-p)(<r)(-«r) 

= 6(-p)(-<r)[(p + cr) - (p)(cr)] - 2(p + tr)[(-p - «r) - (-p)(-<r)] 

= 6(-p)(-a)^(p, a) - 2(p + cr)Si-p, -a) 
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= 6[(-p-cr)-(-p-<r) + (-p)(-«r)]S(p, cr)-2(p + <r)S(-p, -a) 
= 6(-p - <r)5(p, er) - 6/S(-p, -«r)5(p, er) - 2(p + <r)£(-p, -«r) 

-> 4(p + t)S(- P , -a) - 6S(-p, -<r)S(p, a). 
Hence 

I = f° dp (* da- W v (p)W v (o)S(-p, -<r)[4r(p + <r) - QS(p, *)] 

J — oo «'— oo 

which is the form of I used in equation (16) for D (except that P and 
<t are interchanged) . 

The simplification of equations (112) and (116) for T( P , v) and 
U(f), respectively, proceeds along the same lines. In dealing with 
U{f) , the symbolic substitution 

S(p, a, -p, -o-, /) -> [1 + (p)(-p) - 2y' , (-p)] 

•[1 + (*)(-*) - 2y'(-cr)][y r - (/)] 
was found helpful. 

In addition to the simplified forms for T( P , v) and U(f) given in 
equations (17) and (18) , we also have 

S( P) a, v) = S(p + <r, v) - T(p)S(a, v) - T(a)S(p, v). (122) 

XI. THE "SMALL AND SLOW" DEVIATION APPROXIMATION TO We (/) 

This section and the following one are concerned with approxima- 
tions to W s (f) which are obtained by replacing the r's used in equa- 
tions (15) to (20) by the first few terms in their power series expan- 
sions. These expansions are assumed to exist and to converge rapidly 
over the range of frequencies for which W^if) is effectively different 
from zero. Roughly speaking, the top baseband frequency is assumed 
to be small compared with the filter bandwidth. 

When the top baseband frequency is small, the modulating fre- 
quency, /(£), changes slowly and we have the quasistatic case. The 
name "small and slow deviation approximation" is used because (15) 
holds only for "small" rms frequency deviations (D small) , and here 
the requirement of "slowness" is added. 

Two series which play important roles are 

T(/) = 1L«S/n\) a = 1 (123) 

n = 

In T(j) = 2 Kf/nl, X = (124) 
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The first one is the series assumed for r(/). Substituting (123) in 
(124), expanding the logarithm, and equating coefficients of powers 
of / leads to expressions for the A's in terms of the a's: 

X x = a, , 

X 2 = a 2 — a t , 

X 3 = a 3 - 3a 2 a, + 2a\ , (125) 

X 4 = a.i — 4a 3 a! — 3a 2 + 12a 2 ai — 6«i , 

As = « 5 - 5a 4 a! - 10a 3 a 2 + 20a 3 a\ + SOal^ - 60a 2 Oi + 24a* . 

and so on. When the a„'s are the moments of a probability distribu- 
tion, the X„'s are the associated "cumulants" or semi-invariants. In 
our problem the a„'s are proportional to the moments of the normal- 
ized response y(t), a relation which follows when the series (123) for 
r(/) is compared with the one obtained by expanding exp(— i2irjt) in 
the Fourier integral (4) forr(/). 

The small and slow deviation approximation obtained from (15) 
and the first few terms of (123) is 

Wed) -» dl 8(1) 
+ W,(S)[\ + flOiu + 2" 1 Z) 2 X 3 , + 8-Z> 4 X 5 ,) 2 + (a 2r + 2~ l D 2 A r )\] 

+ 2~ 1 (X 2l - + 2"'Z) 2 X 4 ,) 2 f dp W,(p)W 9 (j - p)p\1 - pf 

+ 6" 1 (X a 2 f* dp r da W,(pjWMW ¥ U - p - <r)pV(/ - p- of 

J - J ~ x (126) 

Here the «,,;'s, A,„'s are the imaginary parts of the coefficients in the 

series (123), (124), D 2 = ([/(O/fftr)] 2 ), D is the rms frequency 

deviation in Hz, and a Zr , A r denote the real parts of « 2 , A where 

A = a t — 20:30;! — a\ + 2a 2 o 2 . 

The detailed derivation of equation (126) from equation (15) for 
W g (f) makes use of equation (162) which gives S{x Y , x 2 , • • • x n ) for 
small values of the x's. The leading term in equation (162) gives 

S(p, a) — * pcr\ 2 

S(p, a, v) — » p<n<X 3 

— S(<r, — o-, p, v) + S((x, —a)S(p, v) + 2S(a, p)S( — <r, v) — » a 2 pj/X 4 

(127) 
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where the left side of the last equation is proportional to the integrand 
in equation (112) for T(p, v). There is a similar equation which shows 
that the integrand in the double integral in equation (116) for U(J) 
tends to a quantity proportional to pV/X 5 . To deal with the single 
integral in U(j) we use both terms in equation (162) to obtain 

Sip, -p, /) - ~P%, ~ 2-VfA (128) 

Combining equations (127), (128), and the first three terms in the 
series for T(j) gives 

S(p, a, v) — > pcrvh-i 

T(p, v) -> pv(\ 2 + 2~ 1 D 2 \ i ) 

U(f) -> 1 + fa + 2~ 1 DX + *- l DX)1 + (2"'« 2 + 4T 1 D i A)f. 

Substitution in the small deviation approximation (15) for W g (f) then 
gives the small and slow deviation approximation shown in equation 
(126). 

A form of the small and slow approximation which is more complete 
than (126) may be obtained by starting with the quasistatic form of 
equation (7) for ®(t) instead of from the small deviation approximation 
(15) for W e (t). In the quasistatic case the instantaneous frequency 
q = Wo -f ,p'{t) changes slowly and hence rms cp"(t) is small. This leads 
us to replace <p(t — u) in (7) by the equivalent expression <p(t) — u<p'{t) + 
2 _1 wV'(£) where £ lies between t — u and t. Let F denote the filter 
bandwidth and suppose that the impulse response y(u) is effectively 
outside an interval of length 1/F. Then, heuristically, the integral in (7) 
is given by 

j™ y(u) exp [i<p(t - u)] du = [1 + 0(2~' F~ 2 rms *>")] 

■y(u) exp [i<p(t) — iu<p'(t)] du. 

The integral on the right is the desired quasistatic approximation. It is 
almost equal to the integral on the left when 2~ 1 F~ 2 rms <p" « 1. How- 
ever, for small rms frequency deviations, the contribution of u<p'(t) may 
be less than the term 2 _1 F~ 2 rms <p" even though 

(i) the latter may be <3C 1, and 

(ii) despite the fact that when <p(t) is band-limited with top frequency 
B we always have rms <p" ^ (27rJ5) rms <p'. Therefore, in order to make 
the quasistatic approximation meaningful for small (as well as large) 
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deviations, we impose the additional restriction 2~ 1 F~ 2 rms <p"/(F~ rms 
<p') « 1. Then, if 

rms <p"/(2F 2 ) « 1, rms <p"/(2F rms <p f ) « 1 

the Fourier transform (4) gives the quasistatic approximations 

/ y(u) exp [iip(t - u)] du tt r[</(0/(27r)] exp \i<p(t)] 
J -* (129) 

0(0 « <p(t) - i In r[*/(0/(2ir)] 

which are equivalent to the usual quasistatic approximation for the 
filter output, namely 

s (t) ffrf G(Q) exp [woof + «V(0]. (130) 

D. T. Hess 13 has given a rigorous bound, roughly equivalent to rms 
<p"/(2F 2 ) « 1, for the error in (130). 

For the flat FM baseband case discussed in Section 3.3 the above 
restrictions go into 

10 DB/F 2 « 1, 2B/F « 1 

where D is the rms frequency deviation in Hz and B is the top baseband 
frequency in Hz. Notice that although the term "quasistatic" implies 
that rms <p" tends to in some sense or other, the requirements that the 
quasistatic approximations (129) and (130) hold differ from the require- 
ment that the deviation ratio be large, a condition used in calculating 
the quasistatic approximation to the power spectrum of cos [o) t + <p(t)] 14 
Thus, for the flat baseband case, the deviation ratio can be taken to be 
D/B, and this does not have to be large for (129) and (130) to hold. 
To continue with the derivation of the more complete form of (126), 
we substitute the series (124) for In r(/) in (129) and take the real part. 
This gives 

d(t) « <p(t) + B{i) 
B(t) = E A B ,V(i)/(2*)] B M 

n-l 

Since B(t) depends only on <p'(t), the power spectrum of <p(t) + B(t) is 
W r (j) + W B (j) (Ref. 15). The covariance of B(t) is 

(W -EZ X nf X m ,^' 1 Vr>(27r)- n - m /(n! m!) (131) 

n = l m-1 

where subscripts 1 and 2 denote arguments t and t + r, respectively. 
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From the characteristic function of the joint gaussian distribution 
function of <p[ , <p£ we have 

(gjVT ) _ ,_ „_ 2 f (0+) du [ (0+) <fo exp \-tlu 2 + ^ 2 )2-' - + T w>] 
n\m\ ~ { ' J u J v i n+m u"v m 

[i + (-r-*iri + (-)-*] feM W 



~~ S fc! 4r[(n - fc)2- J + l]r[(m - fc)2" 1 + 1] W„ M2 

(132) 

where ^ T = {<pM, to = (2?r.D) 2 , and we have used, for integer I, 



(friT 1 J 



<o+) „ a' 4- r')(^»/2) 

u-'- 1 exp [-iMY2] <*« = 2r(2-^+l) 



«/a 



Actually, instead of « the upper limit of summation for ft in (132) is the 
smaller of n, ra. Also the sum is unless n, m are both even or both odd. 
When n is even, ft runs over even integers; and when n is odd, k runs over 
odd integers. When (132) is substituted in (131), the Fourier transform 
of the resulting series for (B X B 2 ) gives a series for W B (f) which leads to 
the more complete form of (126) we have been seeking, namely 

■(2t)- 2 " (" # exp (-*2it/t) dr (133) 

The integral in (133) can be expressed as a (ft - l)-fold convolution 
of the power spectrum' IF „.(/) = (2tt/) 2 F „(/) of <p'(t). This gives the 
first few terms of (126), except for the term (a 2r + 2 1 D 2 A r ) W,(j) 
which arises from terms neglected by (133). 

Equation (133) is useful only when D is small because the summation 
with respect to n usually diverges. To illustrate this, consider the single 
pole filter for which T(J) is (1 + ijjTY 1 and X n is (n - 1) K-i/T 1 )"- 
Equation (133), with k replaced by 2ft + 1, gives 

- i ["" (-) n (2w + 2fc)! J D 2n "l 2 

.(2tt)- 4 *- 2 f *" +1 exp (-fiifr) dr. (134) 

The quasistatic approximation for W g (f) obtained by starting with 
(130) -is (see Ref. 13) 
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w.<j) » WM) + g J ^fe r +i)T' /_" *** +1 exp ( ~* 27r/T) dT ' 

/ 24 = / o a 2 * exp[-|--g /,]<&. 

This series, which converges for all D, is of the same form as (134) in 
that it has J 2i /D 2 * +1 in place of the divergent sum with respect to n. 
When D becomes small, the two expressions for We(f) approach equality 
in the sense that the sum with respect to n is the asymptotic expansion 
oiI 2k /D 2h+ \ 

XII. LIOU'S APPROXIMATION FOR SECOND AND THIRD 
ORDER INTERCHANNEL MODULATION 

It is instructive to relate our main result, equation (15) for We(f), 
to an approximation for the interchannel modulation given by Liou. 6 
Liou's approximation is equivalent to taking additional terms in the 
small frequency deviation approximation given in Section XI. 

The interchannel modulation is represented by the portions of 
Wed) in equation (15) which contains T(p, f — p) and S{ P , a, f — p 
— <t). Liou's approximation may be obtained by (i) approximating 
T( p , j — p) by the leading term, namely S( P ,f — p), in equation (17) 
and (ii) expanding S( P , f — P ) , S( P , <t, f — p — <r) in powers of 
p, <r, and / out to and including degree 4. This leads to 

T( P , j - p) m S( P> j - p) = r(/) - T(p)T(i - p) 

= P(f - p)[X 3 + /4 + f% + P (i - P )Q + 0(f) 

S(p, <r, j - p - <t) = pad - p - <r)[A 3 + /A] + 0(f) (135) 

where 

W) = Z*»fM «o = i 

A2 = « 2 — Qfj , X 3 = a 3 — Sc^a, + 2a? 

A = (oii — 2a 3 a x — a 2 + 2a 2 a\)/2 

4 = (a 3 — a 2 ai)/2 

4 = (<*i — <x 3 <*i)/6 

A = 4 _1 (a 4 — a 2 ) — 3 _1 (o;4 — aaaO. 

Equation (162) of Appendix B gives the approximation for S(p, <r, 
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/ — p — a) shown in equation (135). It does not give the higher order 
terms in Sip, f — p) shown in equation (135). These must be cal- 
culated from the series for r(/). Although the <£'s and A's used here 
are not precisely the same as those used by Liou, they are of the 
same general character. 

When the expressions [equation (135)] for T(p, f — p) and S(p, <r, 
f - P - a) are used in equation (15) for Wgif), the second and third 
order interchannel modulation terms are found to be 

I r dp w,(pwm - p)p 2 (i - p) 2 

L J— co 

■{[K i +f%i + P (j- p)t*iT + 1%r\ 

+ I W, + ftl] £ dp £ da 

•WMWMWM - p - *)A 2 (/ - p - *) 2 ( 137 ) 

where the second subscripts r and i denote "real part" and "imaginary 
part." The basic approximation used by Liou [his Eqs. (29) and (30) ] 
may be put in this form by expressing his Fourier transforms as 
convolution integrals and combining terms. 

APPENDIX A 

Power Spectra of Real and Imaginary Parts 
of a Complex Random Process 

Let z (t) be a complex, stationary, ergodic, random process [for ex- 
ample the complex phase angle 0(0] and let x(t), y(t) be its real and 
imaginary parts. We seek convenient expressions for the power spectra 
W x (t), W y (t) of x(t) and y{t) when z(t) is the sum of several correlated 
complex random processes, say a(t), b(t), c(t), • • • . For illustration we 
take 

z(t) = a(t) + b(t) + c(t) (138) 

which corresponds to equation (13) for ®(t) with a(t), b(t), and c(t) 
in place of 3>(t), -iK(t), and i2~ l K 2 (t), respectively. 

Denoting functions with arguments t, t + t by subscripts 1, 2 and 
using relations of the type 

Xl = R C 2l = ( 2l + z *)/2 (139) 

leads to the following expression for the ensemble average <.Tiic 2 )av 
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(X,X 2 )„ = <(Z X + Z*)(z 2 + Z* 2 ))JA 

= <(« A + Z%Z* 2 ) + (z* t Z 2 + 2,2* 2 )) av /4 (140) 

= 2- 1 Re«2 1 2 2 ) av + (^ 2 ) av ). 

It is convenient to write {z x z 2 ) av as 

( 2l z 2 ) av = A(r) + A(-t) (141) 

where 

4(t) = <l(«i«2 + &ib 2 + c^z) 4- Oi&i + aiC 2 4- biC 2 )av . (142) 

This is suggested when the product z x z 2 is multiplied out and terms of 
the type a t b 2 + b^a 2 are considered. Thus, if 

<ai& 2 ) a v = (a(t)b(t + r)> BV = /(r), (143) 

setting t = t' — t and making use of stationarity leads to 

<M 2 )av = (a(t + r)b(t)) av = (a(t'Ml' - r)> av = /(-r) (144) 

and hence to equation (141). 

Similarly, replacing a, , 6, , Cj by a? , 6f , cf leads to writing the 
second ensemble average in equation (140) for (xiX 2 )„ as 

<**z 2 ) av = B(j) + B*{-r) (145) 

where 

B(t) = {\{a*a 2 + 6*6 2 + cfc 2 ) + a*6 2 + a*c 2 + &*c 2 ) BV . (146) 

For terms of the type af b 2 + 6f a 2 , the analogues of equation (143) and 
(144) are 

<a*A>. T = (a*(t)b(t + r)>„ = /(r), 

<6ta 2 ) av = <a(* 4- r)6*(i)>„ = <««')&*(*' - r)>„ 

Comparing equation (13) for ©(£) with equation (138) for z(t) sug- 
gests setting a(t) = *(*), &(<) = ~iK(t), and c(i) = *2~\K: 2 (S); this 
leads to equation (49) for .A(t), B(t) given in Section IV. 

Equation (140) for the autocovariance {x x x 2 )„ of x{t) now takes the 
form 

(x,x 2 ) av = 2" 1 Re [A(t) + 4(-t) + B(t) + B*(-t)], (148) 
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and the power spectrum of x(t) is 

/CO 
exp (-tcor)(.T 1 .r 2 >av dr, co = 2wj. 
-00 

This may be written as 

W £ (f) = 2" 1 Re [P(J) + Q(f) + P*(-/) + Q*(-/)] (1-49) 

where 



P(/) = /* exp (-tor)A(r) dr, 

J-co 

Q(f) = f exp (-tor)B(r) dr. 

J —as 



(150) 



Equation (149) for W x (f) may be derived from 

2" 1 Re [A(t) + A(-t)] - 4-^(7) + 4*(r) + A(-r) + A*(-t)] 
2" 1 Re [B(t) + B*(-r)] = r\B(r) + B*( T ) + fl(-r) + B*(-t)] 

(151) 

and relations of the type 

p*(j) = f exp (wt)A*(t) dr = f exp (-ia>r)4*(-r) dr 
J— 00 «— * 

P(_/) = r exp (-wt)A(-t) dr, (152) 

J-00 

/CO 
exp (-*i>t)A*(t) dr. 
-00 

The power spectrum PP„(/) of the imaginary part y(t) of z(t) may be 
computed in much the same way, starting with 

<»i».>« = «ft - rf)fci - dOWPO" 

= 2" 1 Re [-<8A>„ + <«%>. J. 

This differs from equation (140) for (ziz a )„ only in the sign of <ZiZ 2 )av . 
Therefore only the signs of A(t) and P(/) need be changed in the earlier 
work, and we get 

W v (j) = 2" 1 Re [-P(j) + Q(f) - P*(~f) + Q*(-f)] (153) 
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APPENDIX B 

The Functions S(p, <t),S(p,ct,v), 
The function 

Six, ,x 2 , •■■ ,x n ) = f° duy{u)H u {x x )H u {x 2 ) ■ ■ ■ H u (x n ) (154) 
where 

H u (x) = exp {-i2irxu) - T(x) (155) 

is a symmetrical function of the x's. It may be expressed as the sum 
of products of r's by replacing the H u 's by their definitions, multiply- 
ing out, and using the fact that r(/) is the Fourier transform of y{t). 
This evaluation of the integral shows that S(%, x&, • • • , x n ) is given 
symbolically by 

Six, , x a , • ■ ■ , x n ) = fl [y xt - T(x k )] (156) 

where, after multiplying out, the various powers of y are replaced 
by r's according to the rule y z -> Y (z) . 
For example, 



(157) 



S(p) = 0, 
Sip, a) = [y' - T( P )W ~ r(<r)] 

= y' + ' - y'Tia) - T(p)y' + T(p)T(a) 

= T{p + <r) - T(p)lV) 

Sip, a, v) = Tip + a + v) - Tip + oyI» - T(p + v)T(<r) 

- T(a + v )T(p) + 2T(p)T(<r)T{v). (158) 
For four variables, writing ix) for r(z) , 

Sip, a, v, p) = (p + a + v + p) — (p + a + v)(p) — (p + a + ji)(i») 
- (p + „ + M )(o) - (<r + M + *)(p) + (p + «r)«G0 

+ (p + »)G0G0 + (p + m)WW + (» + ")(p)(m) 

+ (cr + M )(P)W + (" + M)(p)W - 3(p)(o-)(m)W. (159) 

S (a^i , rc 2 , • • • , £„) vanishes when one or more of its arguments are 
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zero because H«(0) is zero. Of interest is the form taken by S(x u 
• • • , x n ) when the x's are small and T (/) may be expanded as a power 
series in /. Let the power series be 



r(/) = £■£«., «o = i. 



n\ 



(160) 



Since 



i2im 



T(/) = f y(u) exp (£/) du, £ = -i 

J —00 

oo en /»oo 

= Zh y(u)?du 

n = "-! J_oo 

it follows that 

« n = / 7(w)r dw. 

•'-00 

When a; is small, equation (155) forff„(aO gives 
H u (x) = exp fee) - T(x) 

= x[({ " «0 + 2" 1 *tf - « 2 )] + 0(x 3 ). 



(161) 



Then 



fl iUO = fox, • • • aO|tt - «O n + 2 -1 (* - «i) n_1 (^ 2 - « 2 ) Z x k 
k-i L i_1 - 1 

+ ooo 

and substitution in the integral [equation (154)] denning S(xy , ■ • ■ , x„) 
leads to 

S(x t , ■■• ,x n ) 
= {x x x 2 • ■ • x n ) J duy(u)\ 2 y^SX-cti)*'' 

+ 2- 1 £ ** % ( n 7 V* 1 " ^ a2)( - ai) """ 1 ] + 0( * n+2) 



— ^XlJ/2 * ' * 3/n/ 






+ 2" 1 E z fc "£ ( W , ^(a,,, - a^X-a,)-'- 1 ] + 0(0. (162) 
This is the approximation needed to examine the form taken by W e (j) 
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when the bandwidth of <p(t) becomes small, as it does in Section XI. 
When T(/) is such that | r(/) - 1 | < e « 1 for all values of /, it 
may be shown that the symbolic form of equation (156) for S(x x , • • • ,x„) 
becomes 

Six, , • • • , x n ) = ft (y Tt - 1) + 0(e 2 ) (163) 

and the T's appear only linearly. Furthermore, S(Xi , • • ■ , x„) is 0(e). 
For example 

S( P , a) = tf + ° - y" - y' + 1 + 0(e 2 ) 

= T(p + c) - r(p)- T(a) + 1 + 0(e 2 ). 

This result is of interest in connection with the first order approximation 
discussed in Appendix D. 

To establish equation (163) let 

z k = y' k - 1, e k = 1 - T(x k ) 
so that equation (156) for S becomes 

S(x> , ■ ■ • , x n ) = fl (** + e,) 

(164) 

= IU+ IK ll'z< + 0(e 2 ). 

k=l 4=1 1=1 

Here the factor z k is omitted from Ip. When the product 

m 

2i«2 ■•■ z m = n (y 14 - !) 

is multiphed out and the y u, s are replaced by r(w)'s, the result is the 
sum of 2 m r's [(1 = r(0)]. Half of the T's will have plus signs and the 
other half will have minus signs. Adding + 1 for each — r and — 1 for 
each +r shows that the entire sum is 0(e). Hence z x z 2 • • • z m is 0(e), 
and when this is used to show that the IP in equation (164) is 0(e), 
the result stated in equation (163) follows. 

APPENDIX C 

The Functions S n and S nm 

The functions S„(xt , x 2 , • • ■ , x„) and «S nm (a; 1 , • • • , x n ; y x , • ■ ■ , y m ) 
are defined by the integrals, for n ^ 1 and m ^ 1-,-. 
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S„(x, , ■■• ,x n ) = / du y{u) exp [o(tt)] II #»(&*)» 

J - M * =1 (165) 

S nm (zi , • • • , x n ; y, , • • • , y m ) 

= / eft; / dwy(v)y(w) exp [a(w) + a(te) + 2c(w, v, o)] 

■'-co »— 00 

Tn ff,(*i)l[ n^)l 

where # u (a;), a(w), c(w, v, o) are given by 
H u {x) = exp (— i2wxu) — T(x) 

a(u) = -\ J" dj W 9 (f)HJL-f)Hj® (166) 

c(w, v, o) = -| f dj W.{j)H w {-f)HXf) = c(v, w, o) 
[see equation (66)]. For n = 0, S is defined as 

S = [ du y(u) exp [o(«)] (167) 

J — oo 

and for the double subscripts, 

SnoOl , • ■ * , %n j) 

/CO /»00 " 

<fo / dw y{v)y(w) exp [o(y) + a(w) + 2c(w, w, o)] II H,(x k ) 
-M •'-00 *=1 

S n0 (xi , ■ • • , x n ;) = £ 0n (; z, , • • • , z„) 

Soo(0 = / ffo / dwy(v)y(w) exp [a(y) + a(w) + 2c(w, v, o)]. 

J-„ J-00 (- 168 ) 

The functions £ n and S nm depend upon both W v (j) and T(j) [through 
H u (f)]. This is in contrast with the function S(xi , • • • , x n ), defined in 
Appendix B, which depends only on r(/) and is independent of W^ij). 

The function S n may be expressed as the sum of multiple integrals 
involving the functions S. Expanding exp [a(u)] in powers of a(u) and 
replacing each a(u) by its integral [equation (166)] with p in place of 
the variable of integration / leads to 

» / INI l"0 fK 

So = 1 + Z Si d * ■ ■ • / dp < 
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r .■ "J 

■ n tt».(p*)js(pi > ■ • • , p, , — pi , • • • , —p>) 

= S(xj , x 2 , • • • , x n ) + 2 ji / dpi ••• dpi 

• n ^(p*) U(pi » • • • i p/ 1 -pi i • • • > -Pi > ^ > • • • > *») 

Lk - 1 J (169) 

where w ^ 1. 

Similarly, expanding exp [2c(w, v, o)] in the integrand of the integrals 
denning the S nm functions leads to 

sm = si+ e ^ r d Pl --- r d Pi 

,=1 ]'• J-co •'-00 

•[n w,Mj8 t (pi , ••■ , Pi)s,(- Pl , ■■■ ,- Pl ) 

&nO\X\ » * " " > #n ]) 

= s s n ( Xl , . • ■ , Xn ) + t, ^ r d Pl ■ . ■ r d Pl . |n w%u) 1 

j = i j ! •'-oo •'-oo La=i J 

•S,+ n (pl , •• • , Pi ,Xl , • • • > -^ n )»S,(— pi , • ■ ■ , — pi) 

S»m(xi , ■ • • , x n ; Vi , • • • , y m ) = Snix,. , ■■■ , x^Smiy! , • • • ,y m ) 

+ E { -=r f d Pi • ■ • P d Pj - [fl tf,(po1 

•S i+n ( Pl , •■■ , Pi ,x x , •■■ , x n )S i+m (- Pl , • • • , -p, , y 1 , • • • , y m ) 

(170) 

when n ^ 1 and m ^ 1. 

In order to obtain an inequality for Snfai , ' ' ' , %n) assume that (i) 
the termwise integration in the derivation of the series in equation 
(169) is legitimate, and (ii) an M > 1 exists such that for all real values 
of/ 

M > | r(/) | = | G(j + /o)/(?(/o) |. 

Since Sfa , • • • , #„) may be expressed as the sum of 2 n terms, each of 
which is a product of not more than n of the r's, 

I S(x, , • • • , x n ) | < TM\ 
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Furthermore, since the integral of W v (/) from — » to + °° is equal to 
(<p 2 (0) B v i the terms of the series in equation (169) for S n (x 1 , ■ • • , x n ) 
are dominated by the terms in 

E ^r (<P 2 )L(2M) 2i+n = (2M) n exp [2MV>.J- 
j=o 3- 

Therefore the series in equation (169) converges and 

I S.(*i ,■■■ ,x n )\< (2M) n exp [2M s (fp a )]. 

This inequahty may be used to show that the series in equation (170) 
for S nm converges and that 

I S na (xi , ■ ■ ■ , *n ; l/i , • • • , y m ) I < (2M)" +m exp [8MVU 

The leading terms in the series required to handle the second and third 
order modulation are 

s = i - \ f dp wms(p, -p) 

+ | /" dp f da W,(p)WM)S(p, a, -p, -a) + 0(/), (171) 
SxO) = - | /" rfp W 9 <j>)S{p, -p, /) 

+ I /" dp F da W *bWMS{p, *, "P, — <r, /) + 0(A (172) 

5 a (p, *) = 5( P| v) - J r rf* W 9 {e)S(*, -<r, p, v) + 0(/), (173) 
S,(p, <r,v) = S(p,a,u) +0&), (174) 

SooCO = 5S + | j^ rfp /" rfcr FT„(p)TF»£(p, ^^-p, -a) 

+ 0(A (175) 

S,o(/;) = &(/) + r rf P P rf * ^,(p)TT,(cr)[4- , S(p, -p, i)S(a, -a) 

J— 00 "—00 

+ 2~ l S(p, f)S(a, -a, -p) + 2~ l S( P , a, j)S(- p, -a)] + 0(<p a ), 

(176) 
&o(p, >0 = S 2 (p, v) - \ f da W p ( ff )S(e, -<r)S(p, v) + Ob'). 
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S 30 (p, tr,v'i) = S 3 (p, a, v) + 0(<p 2 ), 

Sn(p;*) = - f dc WMS(v, p)S(-a, v) + 0(<p 4 ), 

J —DO 

S 2 i(p, <r',v) = 0(<p 2 ). 

In obtaining the leading terms in S o(',) and $ ]0 (/;) the leading terms in 
So , $i(/) and S 2 (p, /) were used. *S t (/) is 0(<p 2 ) in contrast with S n , 

S 2 , S 3 , ■ • • which are 0(1). 

APPENDIX D 

Derivation of Earlier First Order Approximation 
by Present Procedure 

The first order approximations which are given in Section VII 
are somewhat more complicated, as well as more accurate, than the 
ones which have appeared in the literature. 3, *■ 0_12 Here the relation 
between the earlier and present work will be brought out by applying 
the procedure of Section VII to obtain a first order approximation, 
which is the same as the given in Ref. 10 [the 6 it) of Ref. 10 is 9{t) 
— f(t) in the present notation]. 

The derivation starts from the initial equations (5) and (7) for 
the filter output So(t) , 

s (t) = {exp [— a — ip + i®(t)]\ exp (ico t) . . 

©(0 = — i In < / du y(u) exp [i<p(t — u)] f • 

The difference between the output phase angle 8(f) = Re S(t) and the 
input phase angle <p(t) is assumed to be small, and the filter delay is 
usually taken to be zero at the carrier frequency, that is, Im [dT(f)/df] 
is zero at / = 0. 

Adding and subtracting i<p(t) [instead of the linear portion $(t) of 
the output] in the exponent appearing in equation (177) for @(t) gives 

@(t) =<p(t) -tln[l +*«], (17g) 



k(t) = I duy(u)\exp [i<p(t — u) — i<p(t)] 

v — er\ 



11. 



The first order approximations to the complex phase angle ©(0 and 
the output phase angle 6(t) are now 

0(0 pa <p{t) - ik(t) 
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and 

0(0 « <p(t) + Re [-ik(t)], (179) 

respectively. 

The analysis of the earlier sections goes through much as before with 
k(t) in place of K(t), and 

/*„(/) = exp (-i'lTJu) - 1 (180) 

in place of 

ff-CO = exp (-fiwfu) - r(/). 

An illustration of how fc u (/) enters the analysis is furnished by the 
computation of (exp [vp(t — u) — i<p(l)]) av . As suggested by equation 
(56) for <p(t - u) - * (0,1st 

L(p(t) — <p(t — u) — <p(t) 

€(f) = exp (-12T/U) - 1 = />„(/) ( 

-| /_* df W.Wm-f) = ~\ £ df W v (j)h(f)K(-f) = a'(u) 

(exp \itp(t — u) — «V(0]>nv = exp [a' (it)]. 

Of most interest in practice is the power spectrum W^(j) of £(0, 
£(0 = Re [-ik(t)} (lg2) 

e(0 « <p(t) + £(0 

where £(0 is an approximation to the distortion. The power spectrum 
Wt(j) is the Fourier transform of the covariance (£i£ 2 ) where, as before, 
subscripts 1, 2 refer to times t, t + t, respectively. By putting £(0, 
— ik(t) for 0(0, ©(0 in equation (47), or directly, 

&&>« - 2" 1 Re [-<*A>« + <*#.>„]. (183) 

It may be shown that 

(fci)av = / du y(u) exp [a'(u)] 

■» -oo 

<fcifc 2 )av = <fci) n v + ( dM f <fo t(w)7(«) exp [a'(u) + o'(»)] 

•{exp[2c'(u,», r)] - 1} 
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(k*M» = | (fci)av | 2 + [ du f dvy*(u)y(v) exp [a'(u) + a'(v)] 

•'-00 •'-00 

•|exp[-2c'(u,y, t)] - 1} 
(184) 
where 

c'(u, v, r) = -\ r df W,(j)h u (-f)h v (f) exp (;2tt/t) 

a'(u) = c'(«, m, o) (185) 

c'(u,v, —t) = c'(v,u, t). 

Since h(—f) is equal to /&*(/), a'(w) and c'(w, v, t) are real. Furthermore, 

c'(u, v, t) = -$[fl„(r + u - v) - R^t - v) - R v (t + u) + #,(r)] 

a'(u) = -R v {o) + fl» (186) 

where R v (t) is the covariance (<p(t)<p(t + t)) ov of «?(0- 

The expression for (£i£ 2 )ar obtained by combining equations (183) 

and (184) is similar to equation (8) of Ref. 10. 
The power spectrum of £(£) may be written as 

W t (J) = & «(/) + 2" 1 Re [P(f) + Q(f) + P*(-1) + Q*(-f)] (187) 

in which £,„ is equal to Im <A;,) av and 

P(/) = P dr exp [-i2 7 r/r][-K(fc 1 fc 2 ) BV - (fc,)^)], 

J -° (188) 

Q(/) = r dr exp (-;27r/r)[K<Wav - I (fci)av | 2 )]. 

J_oo 

Addition gives 

/oo /»oo /•« 

rfr exp (-i2tt/t) / dw / dv exp [a'(u) + a'(v)] 
-00 J — 00 «'—00 

■(-h(u)y(v){exp [2c'(u,v, r)] - 1| 

+ §Y*(«)y(») I exp [-2c'(u, v, t)] - 1}) (189) 

which, when used in equation (187), leads to an expression for W^(j) 
which is similar to the main result given in Ref. 10 [equation (16) of 
Ref. 10]. 
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The relation c'(u, v, — r) = c'(v, u, t) may be used to show that Q(f) 
is real and P(/) is even. When y(u) is real, T(— /) = T*(j), and the 
expression for W^ij) may be simplified. 

Expanding exp [±2c'(w, v, t)] in powers of c'(u, v, t) leads to 

Wt(f) = & 5(/) + T S ^ T dP. • • • f ^n «(/ - Pi - • • • - P.) 

i n=l '*' "-» «'—oo 

f[ W%U) | «(p, , • • • , P») - (-)'«*(-* , • ■ • , -P„) I 2 

k=l 

(190) 
where, as in Appendix C for the unprimed S's, 
S&i , • • • ,*n) 

/.oo n 

= / duy(u) exp [a'(w)] H t ex P ( — *2irtia;*) — 1] 

J -00 *=1 

= S'fo , • • • , aO + Z m } - r d Pl • • • P dp, |n WJjfi 

* S'(pi , • • • , (>j , — Pi , • • • , — Pj , X\ , • • • , Z„) . 

The series in equation (190) is analogous to the series in equation (90) 
for the more accurate first order approximation based on <p(t — u) — 
$(t). The function S'fa , • • • , x n ) is the analogue of S(xi , • • • , x n ) 
discussed in Appendix B and is defined by 

S'(*. , ••■ ,x n ) = r duy{u) f[h u (x k ) 

•'-oo i = l 

= nor'-i). 

A = l 

The second equation is symbolic in that ?/* is to be replaced by T(z) 
after expansion of the product. It is shown in Appendix B that when 
| T(/) - 1 | < e « 1 for all real values of /, 

S(x t> ■■■ ,.T„) = S'(x 1} ■■■ ,x n ) +0(e 2 ). 
APPENDIX E 

Power Spectrum of ad(t) + b In R (t) 

Section III states that an expression for the power spectrum W x (f) 
of x(t) = ad(t) + b In R(t) may be obtained from equation (15) for 
Wed) by replacing U(f), T(p, f-p), and S(p, <r, j-p-tr) by (a+ib)U(f), 
(a + ib) T(p, / — p), and (a + ib)S(p, <r,j — p — a), respectively. Here 
the steps leading to this result are outlined. 
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From equation (6) for the complex phase angle ®(i) we have 

i8(0 = In R(t) + i@(t), 

and it follows that, for arbitrary real values of a and b, 

ad + b In R(t) = Re [(a + ib)@(t)]. 

Consequently W x (j) may be obtained by replacing 0(0 by (a + ib)0(t) 
in the analysis which led to equation (15) for W B (f). 

The functions A(t) and B(t) appearing in equation (48) are replaced 
by (a + ib) 2 A(T) and \a + ib\ 2 B(t), and their respective Fourier 
transforms P(f) and Q(f) are replaced by (a+z'6) 2 P(/) and | a+ib | 2 Q(f). 
Each factor in (a + ib) 2 = (a + ib)(a + ib) can be associated with 
factors in P(j), and each factor in (a + id) (a + ib)* with factors in 
Q(j), in such a way that £/(/) becomes multiplied by (a + ib) and 
U*(—j) by (a + z7>)*, and so on. This may be verified by repeating the 
analysis of Sections VIII and IX with the modified expressions. 

APPENDIX F 

Results Obtained from the Series for ln[l + K(t) ] may be Asymptotic 

For gaussian <p(t) with average and rms value <r, the following con- 
siderations suggest that results obtained from equations (10) and (12), 
namely 

©(0 = $(t) - i In [1 + K(t)] (191) 

= $(0 + i T d n l [-K{t)] n 
i 

represent the first few terms of an asymptotic series when a — > 0. 

Since K(t) is difficult to handle, we replace it by a[a l,eU) — 1] where a 

is somewhat like the integral of y(u) between — » and oo. The value 

of this integral is 1, and we regard a as being near 1. The series for 

In [1 + K(t)] behaves somewhat like the series 

In (1 + a{exp Mt)] - 1}) = - Z ( -~^ fexp [vp(fi] - l} n (192) 

i n 

in which the mean square value of the modulus of the nth term is 

f 6XP [ ~ y Y/ 2<r2 )] n~\2a sin ^/2] 2 " d<p. (193) 

When o- <3C 1 and n is not too large, most of the contribution to the 
integral (193) arises from the region around <p = 0, and the integral is 
approximately 

1-3 ■■■ (2n- l)(aa) 2n /n 2 . 
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Consequently, the first few terms decrease rapidly when <r <K 1. 
However, when n is very large most of the contribution arises from the 
regions around p = ±ir, ±3ir, • ■ • , where [sin (<p/2)] 2n is a narrow pulse 
of height 1 and area 2(ir/n) i . When a « 1 only the regions around 
<p = ±7r are important and the integral (193) is approximately 

(r -l 2n 3/2-5/2 (2a)2 „ exp [ _ 1r 7( 2(r 2 )] 

which tends to °o when a is near 1 and n — » » . 

The fact that the rms values of the terms of the series in equation 
(192) decrease rapidly at first and then increase without limit suggests 
that results attained from the somewhat similar series in equation (191) 
may be asymptotic in nature as a — * 0. 
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